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Abstract 



The collective modes of QCD at temperatures and densities above its phase- 
transition are analyzed; models for various observables based on the knowledge 
gained from this analysis are constructed, and the results derived are compared 
to those of other methods, where available. Specifically, isotropic systems as well 
as systems having an anisotropy in momentum-space are investigated using the 
Hard-Thermal-Loop (HTL) approximation. Several observables are calculated, 
ranging from thermodynamic quantities and their application in the isotropic 
case to the collisional energy loss in the isotropic as well as the anisotropic case. 

For isotropic systems, results at finite densities (or non- vanishing quark chem- 
ical potential) based on a phenomenological description of lattice calculations at 
vanishing densities are derived and compared to those following from different 
methods. It is shown that results for the thermodynamic pressure and entropy 
as well as those for the quark-number susceptibilities agree well with indepen- 
dent lattice calculations, and it is demonstrated that the plasmon effect leads 
to numerically small contributions, in contrast to what is found in strictly per- 
turbative approaches. As possible applications of the resulting equation of state 
the mass-radius relationship of so-called quark-stars as well as the hydrodynamic 
expansion in the Bjorken model of the quark-gluon plasma created through a 
heavy-ion collision are calculated. 

It is shown that systems with anisotropic momentum-space distributions con- 
tain unstable modes in addition to the stable quasiparticle modes, which - since 
they correspond to exponentially growing field amplitudes - may be of great 
importance for the dynamical evolution of an incompletely thermalized quark- 
gluon plasma. Moreover, the presence of such instabilities and the correspond- 
ing singularities in the propagator lead to divergences of scattering amplitudes 
in a perturbative framework, signalling the breakdown of the latter. However, 
it is demonstrated that at least one observable, namely the collisional energy 
loss, is protected from these divergencies. This permits its calculation also for 
anisotropic systems, as can be shown analytically both for very weak as well 
as for extremely strong anisotropies. A subsequent comparison of results from 
isotropic and anisotropic systems exhibits a possibly strong directional depen- 
dence of the energy loss for the latter, which might lead to effects that can be 
verified experimentally. 
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There is a theory which states that if ever anyone 
discovers exactly what the Universe is for and why it is 
it will instantly disappear and be replaced by 
something even more bizarre and inexplicable. 



There is another theory which states 
that this has already happened. 

Douglas Adams, The Hitch Hiker's Guide to the Galaxy II 



Chapter 1 
Introduction 



1.1 Quantum chromo dynamics 

Today, quantum chromodynamics (QCD) is accepted to be the established theory 
of strong interactions. It has been formulated along the lines of the theory of 
quantum electrodynamics (QED), which - as the unification of quantum theory 
and electrodynamics - is certainly one of the most successful and accurate theories 
in modern physics. Consequently, QCD bears several similarities but also exhibits 
some striking differences with respect to QED, as will be briefly illustrated in 
the following. Both theories are gauge field theories, but whereas QED is an 
Abelian gauge theory with the gauge group U(l), quantum chromodynamics is 
non-Abelian in nature, having the color group SU(3) as a gauge group. The 
gauge boson associated with the U(l) group of QED is the well-known photon, 
while for the color group SU(3) there are 8 associated gauge bosons called gluons. 
In contrast to the photon which is uncharged, the gluons do carry color charges 
which are the QCD equivalents of the electromagnetic charge. For QED, the 
matter particles are the electron, muon and tauon while for QCD there are the 
six quark species up, down, strange, charm, beauty and top, which are all spin 
1/2 fields or fermions; those quark species (or flavors) with masses much smaller 
than the energy scale under consideration are referred to as active or light flavors 
and their number is usually denoted by Nf. 

Both QED and QCD are renormalizable field theories and as a consequence 
their coupling "constants" are functions of the energy scale Q. In Abelian gauge 
theories (as QED) the coupling increases for larger energy scales, whereas for low 
energies it turns out to be small, allowing a valid (and indeed extremely accurate) 
description of our (low energy) everyday world using perturbation theory. 

However, unlike in Abelian gauge theories, the non-Abelian nature of QCD 
makes its coupling constant a s (Q) decrease as Q becomes large, a property that 
is known as asymptotic freedom. Consequently, at very high energies one may 
hope that perturbation theory offers a valid description of QCD since the coupling 
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constant gets tiny. On the other hand, when the energy scale gets smaller and 
the coupling constant rises, calculating observables in QCD becomes very hard 
in general, since one has to deal with a strongly coupled theory and perturbation 
theory breaks down. Therefore, one obviously needs non-perturbative methods 
to describe the theory and explain the experimental fact that at very low energies 
quarks and gluons cannot roam freely but are confined into hadrons from which 
they do not escape. 

While at low temperatures quarks and gluons are thus locked up mainly inside 
protons and neutrons, one expects them to propagate freely in a state that has 
been dubbed the quark-gluon plasma at very high temperatures. The exact nature 
of this transition between the confined phase and the quark-gluon plasma phase 
(whether it is a first or second order phase transition or a crossover) is still a 
matter of active research. 

Several theoretical methods have been proposed to study the physics of QCD, 
each most apt for a certain energy range, the most popular being: 

• Effective models based on the QCD Lagrangian that allow to calculate the 
low-energy behavior of the theory. Perhaps the most interesting of these 
is chiral perturbation theory, based on the effective chiral Lagrangian (see 
Gasser's and Leutwyler's original articles on this subject [1, 2] or [3] 
for a pedagogical review). 

• Monte-Carlo simulations of QCD on a lattice, which so far is the only 
method to provide quantitative results for the intermediate energy range 
near the deconfinement transition. Since simulations with dynamical quarks 
are much more time consuming than those without, the situation best stud- 
ied in lattice QCD is that of a pure gauge theory. Nevertheless, there has 
also been considerable progress for systems with quarks and recently even 
for non-vanishing chemical potential. The latter in general poses a major 
problem for lattice calculations because of the sign problem, which pro- 
hibits the use of conventional numerical algorithms (see Karsch [4] for a 
pedagogical review). 

• Perturbative approaches, which are expected to be most accurate for very 
high energies. For instance, there exist several methods to calculate the 
equation of state for the quark-gluon plasma phase; a short review on this 
subject will be given in section 1.4. 

On the experimental side, the main earth-based tool to investigate QCD near 
or above the phase transition is ultrarelativistic particle collisions. The study 
of these makes huge collider facilities necessary, such as those of CERN {Cen- 
tre Europeenne pour la Recherche Nucleaire) in Switzerland/France and BNL 
(Brookhaven National Laboratory) in the U.S.A. While the SPS (Super Proton 
Synchrotron) collider at CERN produced the first experimental indication of the 
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existence of the quark-gluon plasma phase, it has currently been dismantled to 
make way for the next generation collider LHC (Large Hadron Collider), which 
will allow to probe energies well above the phase transition and is due to be- 
come operational in 2007. Since 1999, RHIC (Relativistic Heavy Ion Collider) 
is the collider operational at BNL, which has already produced a wealth of new 
data bringing the study of nucleus-nucleus collisions carried out at SPS to a new 
energy regime. 

Finally, indirect studies of QCD may be possible through astrophysical obser- 
vations and cosmology [5]. 

1.2 The deconfinement transition 

Our current understanding of QCD is such that whenever the system under con- 
sideration becomes sufficiently hot and/or dense, hadrons dissolve into a gas of 
almost free quarks and gluons, the quark-gluon plasma. More precisely, either 
the temperature T, or the quark chemical potential /x, or a combination of both, 
has to be big enough so that the system undergoes a transition to the deconfined 
state of matter. In a slight abuse of terminology, one refers to the temperature 
and chemical potential, for which the system undergoes the transition, as critical 
temperature T c and critical chemical potential /i c , regardless of the exact nature 
of the transition (first/second order phase transition or crossover). 

At vanishing chemical potential, lattice studies of QCD tell us that the critical 
temperature is on the order of a few hundred MeV, which is much smaller than 
the masses of the charm, beauty and top quark. For this reason, those heavy 
quarks do not play a role in the description of the physics near the deconfinement 
transition, so it is useful to limit the number of active (=light) quark flavors to 
Nf = 3. Moreover, since the up and down quark masses are so small compared 
to the energy scale of the deconfinement transition, these quarks may be taken 
to be massless. The strange quark, finally, is neither much heavier nor much 
lighter than the deconfinement scale, so the number of active quark flavors near 
the deconfinement transition is between two and three, also sometimes denoted 
as Nf — 2 + 1. 

Based on universality arguments, it is possible to show that for a pure gauge 
theory (Nf = 0, corresponding to infinitely heavy up, down and strange quarks) 
as well as for three and more active flavors (Nf > 3, corresponding to zero 
up, down and strange quark masses) the deconfinement transition should be a 
first-order phase transition. For Nf = 2, the transition is probably a second 
order phase transition, but the situation is not entirely clear [6]. Fig. 1.1 shows a 
cartoon of our current knowledge about the nature of the deconfinement transition 
at vanishing chemical potential, indicating that for light up and down quarks 
and a moderately heavy strange quark the transition is probably a cross-over [7], 
though a weak first-order transition cannot be ruled out either (if the physical 
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pure gauge: 
2(3} 






chiral limit: 
sum x SUiZ) 



Figure 1.1: Character of the transition as a function of the light quark masses 
m u = ma and the strange quark mass m s . The dot represents the physical values 
of the parameters. From [6]. 



point in Fig 1.1 moves down). 

The value of the critical temperature has been studied intensively in lattice 
QCD calculations at \i = 0, giving [8, 4] 

• T c = 271 ± 2 MeV for N f = 

• T c = 171 ±4 MeV and T c = 173 ±8 MeV for N f = 2, depending on whether 
one uses Wilson or staggered fermion lattice implementations 

• T c = 154 ± 8 MeV for N f = 3. 



1.3 The quark-gluon plasma 

The quark-gluon plasma phase is unlike any of the other phases we know of in 
nature: it is a new state of matter, on equal footing with solid, fluid, gaseous and 
the electromagnetic plasma, which makes its study interesting from a conceptual 
point of view already. Furthermore, from standard cosmology we expect that the 
quark-gluon plasma was the prevailing form of matter until about 10~ 5 seconds 
after the Big Bang, which could have interesting implications on the universe. 
Nowadays, it is expected to be produced only in ultrarelativistic heavy-ion col- 
lisions, and perhaps reside in the core of heavy neutron stars, which might have 
interesting astrophysical consequences. For these reasons, one would like to know 
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more about the properties of this state of matter, especially since ultrarelativistic 
heavy-ion collision measurements in the next decade may allow to verify these 
predictions experimentally. 

One of these properties is the equation of state (EOS) of the quark-gluon 
plasma; once the latter is created through the collision of two ultrarelativistic 
nuclei, the expansion of the "fireball" is controlled mainly by the EOS. Also, 
the EOS specifies the mass-radius relation of so-called quark-stars through the 
Tolman-Oppenheimer-Volkoff equations. 

The EOS itself follows from the pressure p(T, jj) (which is related to the ther- 
modynamic potential Vt = —pV) through the standard thermodynamic equations 

dp dp 

s = -, n = -, e=-p + sT + n^ (1.1) 

where s, n and e are the entropy, number and energy density, respectively; from 
these, the energy density is given as a function of the pressure only, e = e(p), 
which is one representation of the EOS. As a consequence, knowledge of the 
pressure for arbitrary temperature and chemical potential automatically implies 
the EOS. 



1.4 The QCD pressure 

Since the QCD coupling becomes small because of asymptotic freedom, one could 
expect that weak coupling calculations at high temperature T or chemical poten- 
tial \i lead to reasonable estimates of the thermodynamic pressure of QCD. To 
zeroth order in the coupling (where the quark-gluon plasma would correspond to 
an ideal gas of quarks and gluons), the pressure takes the form 

(A 2 -l)vr 2 T 4 ArAr /7vr 2 T 4 /i 2 T 2 // \ 

where Nf is the number of massless quark-flavors and A = 3 is the number of 
colors [9]. The first correction to this result (which itself corresponds to black- 
body radiation in electrodynamics) is of first order in a s and given by 



36 v " ~'~ 1 144 V"~ ' ~~ ^ 2 ff2 

(1.3) 

1.4.1 Vanishing chemical potential 

To simplify the argument, consider the chemical potential to be zero for the 
moment, so that the pressure is proportional to T 4 and all the non-trivial tem- 
perature behavior of p/T 4 resides in the temperature dependence of the strong 
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coupling a s . The standard procedure is to take the solution for the running cou- 
pling from the renormalization group equation, which to lowest (1-loop) order is 
given by 

12tt 

asW ~ (UN-2N f )lnfiyAy {1A) 

and put the renormalization point p (not to be confounded with the chemical 
potential /z) proportional to the first Matsubara frequency, p = c^nT, which for 
massless quarks is the only dimensionfull quantity inherent to the theory. The 
parameter expresses our ignorance about the ideal renormalization point and is 
usually varied by a factor of two to obtain an estimate of how strongly the result 
depends on this chosen renormalization point; obviously, one would prefer the 
result to have as little dependence as possible on the arbitrary parameter c^. The 
strong coupling a s also depends on the scale parameter of the modified minimal 
subtraction scheme, A^g, which cannot be fixed by theoretical considerations 
but has to come from experiments. However, a s is not easily measured at energy 
scales near the deconfinement transition, so one has to extrapolate measurements 
down from higher energies. Currently, the official extrapolation value at p — 2 
GeV is given by [10] 

a s (p = 2GeV) = 0.2994; (1.5) 

since from the choice of p ~ 2nT this corresponds to a temperature where one 
expects up, down and strange quarks to be effectively massless (Nf = 3), this 

(3) 

would fix Aj-^ ~ 194MeV. However, by using a higher order correction to the 
running coupling (2-loop), 

„ (rA = 12n fi - (^-13NN f + 3N f /N)\nL(p) \ 

sW (UN-2N f )L(p)\ 6(llN-2N f yl(p) J [ ' 

where L(p) = In p 2 / A^^, the value of A^jg changes to A^- ~ 378MeV and adopt- 

(3) 

ing a 3-loop running coupling, one finds A^- ~ 344 MeV. One could still improve 
on this result by taking into account matching conditions between the different 
quark mass scales [11] as well higher order (e.g. 4-loop [12]) beta-functions for 
the running coupling; however, being interested in a rough estimate only, one 
can naively use the critical temperature T c ~ 154MeV for Nf = 3 to obtain the 



ratio T c / A^-L ~ 0.45, which turns out to be reasonably close to the value for two 



MS 

active flavors advocated by Gupta [13] 

T, 



^0.49. (1.7) 
MS 



Therefore, I will adopt this ratio for all the two-flavor calculations in the 
following. Clearly, if a different value is adopted, some quantitative results in 
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1.5 2 2.5 3 3.5 4 

T/T c 



Figure 1.2: Strictly perturbative results for the thermal pressure of Nf = 2 QCD, 
normalized to the ideal-gas value p , as a function of T/T c (assuming T c /A^jg = 
0.49). The various gray bands bounded by differently dashed lines show the 
perturbative results to order a s , a?J 2 , o? s and a J 2 with the MS renormalization 
point fl varied between nT and 4nT [17]. The boxes and triangles show the 
estimated continuum-extrapolated lattice results from [18, 19]. 



this work will change somewhat, but the qualitative picture remains the same 
and the modifications to the results should be straightforward (I will also give 
the trend of the modified result in chapter 5). 

Having determined the temperature dependence of the coupling one can pro- 
ceed to evaluate the perturbative pressure; indeed, explicit calculations of higher 
order contributions to the pressure have been pushed to aj 2 [14, 15] and even to 
a\\o^a s in a recent heroic effort [16]. The result, however, is rather depressing: 
instead of showing any sign of convergence, adding successive orders gives a total 
pressure that is sometimes below, sometimes above the free pressure po, jumping 
in an nearly unpredictable fashion (see Fig. 1.2). Perhaps worse, the ambiguity 
introduced by the constant seems to increase for every order added to the 
pressure, signalling a complete loss of predictive power. 

Various mathematical extrapolation techniques have been tried, such as Pade 
approximants and Borel resummation [20, 21, 22, 23, 24]. The resulting expres- 
sions are indeed better behaved than polynomial approximations truncated at 
order aj 2 or lower, showing a weaker dependence on the renormalization scale. 
However, while these methods do improve the situation somewhat, it is fair to 
say that they offer little physical insight on the source of the difficulty. 

However, in recent years resummations based on the hard-thermal-loop (HTL) 
effective action [25] have been proposed, alternatively in the form of so-called HTL 
perturbation theory [26, 27] or based upon the 2-loop ^-derivable approximation 
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[28, 29, 30] (see also Peshier [31]) 1 . The latter approach, which assumes weakly 
interacting quasiparticles as determined by the HTL propagators and NLO cor- 
rections thereof, leads to results which agree remarkably well with lattice data 2 ; 
more details about this approach will be given later on. Recently, it has also been 
shown that the results from the HTL quasiparticle models are in fact consistent 
with those obtained in high-order perturbation theory if the latter is organized 
through an effective (dimensionally reduced) field theory according to [35] and 
effective-field-theory parameters are kept without further expanding in powers of 
the coupling [17], as already advocated in Ref. [16]. 

1.4.2 Finite chemical potential 

In contrast to vanishing chemical potential, where lattice studies seem to be 
quite capable of calculating the EOS for a quark-gluon plasma near the decon- 
finement transition, this task is by far more difficult at finite (and especially large) 
chemical potential. The reason for this is that the quark determinant becomes 
complex at finite chemical potential /x, so Monte-Carlo simulations cease to be 
directly applicable [36]; this is known as the sign problem. Recently however, 
there has been important progress also for lattice calculations 3 at non-vanishing 
/i [38, 39, 40, 36, 41, 42, 43], which circumvent the problem of a complex quark 
determinant by either a Taylor expansion in \i of the observables at vanishing 
chemical potential or an analytic continuation of observables calculated at purely 
imaginary \i. Unfortunately, all these calculations are limited to the case of small 
chemical potential (fJ>/T <C 1) for now, so an EOS for cold dense matter, which is 
of importance in astrophysical situations [44, 45, 46, 47]) is currently beyond the 
reach of lattice calculations. Moreover, a recent result for a model of QCD with 
an infinite number of flavors (Nf — > oo) indicates that extrapolation of results 
obtained at \i = breaks down rather abruptly for fi/T > ir, pointing presumably 
towards a generic obstruction for extrapolating data from small to large chemical 
potential [48]. 

As a remedy for this situation, Peshier et al. [44, 45] proposed a method 
which can be used to map the available lattice data for /j, = to finite \i and 
small temperatures by describing the interacting plasma as a system of massive 
quasiparticles. Since this procedure does not involve any direct extrapolation, the 
model proposed by Pehsier and its hard-thermal-loop (HTL) extension which 

1 For recent progress with higher-loop ^-derivable approximations in scalar field theory see 
Refs. [32, 33, 34]. 

2 This is in fact not the case for 2-loop HTL perturbation theory [27]. However a dimension- 
ally reduced variant [17] of this form of improved perturbation theory (which is perturbatively 
equivalent and moreover much simpler) has recently been shown to agree well with the results 
of Refs. [28, 29, 30] and correspondingly also with lattice data. 

3 A perturbative expansion of the pressure for small chemical potential has also been calcu- 
lated recently [37]. 
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will be considered later on are generally capable of calculating the EOS for the 
quark-gluon plasma for any (i.e. also large) value of fi/T, given that the models 
stay reliable for all temperatures and densities. Accordingly, such quasiparticle 
models 4 are a valuable complementary method to e.g. lattice calculations at 
small chemical potential, and they turn out to give quite robust quantitative 
predictions on the EOS for cold dense QCD matter as well as on the critical 
chemical potential. 

I will give a detailed discussion of the results for the EOS obtained by these 
quasiparticle models and also compare with results from other approaches [51, 
52, 46, 47]. 

1.5 Anisotropic systems 

So far, all the discussion has been limited to homogeneous and isotropic systems. 
However, considering that the system created through a collision of two heavy 
ions propagating essentially at the speed of light is certainly not homogeneous 
and isotropic, one has to ask whether a description through equilibrium field 
theory actually constitutes a valid treatment of this system. This question can 
be restated in the form of whether the system thermalizes fast enough so that at 
least at a local scale it can be described by equilibrium field theory. The advent 
of "bottom-up" thermalization by Baier, Muller, Schiff and Son [53] was a 
big step towards answering this question, although there are indications that this 
scheme probably will have to be modified because of more recent developments 
discussed below. 

Despite continuous progress [54, 55, 56, 57], an out-of equilibrium study of 
QCD has not been achieved yet; however, one may try to learn more about 
out-of-equilibrium situations by studying systems that are close to equilibrium. 
For instance, one can consider the effects that an anisotropic momentum-space 
distribution function has on the thermalization, or more precisely isotropization 
of a quark-gluon plasma. As has been pointed out first by Mrowczynski [58, 59, 
60], such a system possesses instabilities that correspond to so-called Weibel or 
filamentation instabilities in electrodynamics [61]. Weibel showed that unstable 
transverse modes exist in electrodynamic plasmas with anisotropic momentum 
distributions and derived their growth rate in linear response theory. In plasma 
physics, it has been shown that these instabilities generate strong magnetic fields 
resulting in the filamentation of the electron current by simulations and recently 
also experimentally [62]. The effects of these instabilities are much less clear, but 
they may potentially be very important for the quark-gluon plasma evolution at 
RHIC or LHC due to the large amount of momentum-space anisotropy in the 
gluon distribution functions at r ~ 1 fm/c. 

4 See also the model in Ref. [49] which has been extended to finite chemical potential in 
Ref. [50]. 
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Mrowczynski and Randrup recently performed phenomenological esti- 
mates of the growth rate of the instabilities for two types of anisotropic dis- 
tribution functions [63], finding that the degree of amplification of the Weibel 
instability is not expected to dominate the dynamics of a QGP, but instead is 
comparable to the contribution from elastic Boltzmann collisions. However, if a 
large number of the unstable modes could be excited then it is possible that their 
combined effect on the overall dynamics could be significant. 

Moreover, Arnold, Lenaghan and Moore [64] argued that the presence 
of instabilities giving rise to large gauge fields necessarily change the first stage 
of the bottom-up thermalization scenario qualitatively, making a revision of the 
original scheme unavoidable. 

In this work, I will perform a detailed study of the HTL quasiparticles and 
unstable modes in an anisotropic quark-gluon plasma and discuss their effects on 
the properties of such a system with respect to an isotropic quark-gluon plasma. 

1.6 Outline of this work 

In chapter 2, which is still of introductory character, I will give a definition of 
what I will refer to as quasiparticles and also introduce a simple quasiparticle 
model for the QCD pressure of the isotropic quark-gluon plasma. In chapter 3, 
I will introduce resummed perturbation theory based on the 2PI effective action 
and an extension of the simple quasiparticle model that follows from it. In chapter 
4, I will discuss the results one obtains from mapping lattice data from zero to 
arbitrary chemical potential potential, including two possible applications. In 
chapter 5, I will investigate to what extent these results are modified when one 
limits the input to only the value of T c /Ay^. In chapter 6, I will calculate the 
HTL resummed gluon self-energy of an anisotropic quark-gluon plasma, identify 
all stable and unstable collective modes of the system, and consider their effect on 
the partonic energy loss with respect to an isotropic system in chapter 7. Finally, 
I will give my conclusions in chapter 8. 



Chapter 2 

QCD Quasiparticles 



When elementary particles propagate through a medium, their vacuum properties 
are changed through the effect of their interactions, which is also referred to as 
"dressing". In the simplest case, a quasiparticle is then just an elementary particle 
which had e.g. its vacuum mass changed to an effective mass by the medium. 

However, in general a medium is characterized by a whole set of collective 
modes or quasiparticles (I will use the terms interchangeably), some of which cor- 
respond to elementary particles in the vacuum and others which do not. Loosely 
speaking, a quasiparticle is an in-medium excitation that behaves as if it was an 
elementary particle, although it might not correspond to one in the vacuum. 

Generally, quasiparticles possess a dispersion law E(k) that gives their energy 
E as a function of their momentum k and - since their lifetime may not be infinite 
- a decay or damping rate j(k). Both the dispersion relation as well as the decay 
rate are linked to the peak and width of the spectral densities; in the simple 
cases considered in the following the spectral densities can be assumed to have a 
Breit-Wigner form so that the dispersion relations and decay rate are linked to a 
simple pole of the propagator. 



2.1 The gluon propagator in the HTL approxima- 
tion 

The gluon propagator G^ v is related to the gluon self-energy U^ u through the 
Dyson-Schwinger equation 



gJ, — g^ + n 



[IV 



11 



(2.1) 
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Using a tensor decomposition of the self-energy and restricting to the HTL ap- 
proximation [65, 29] and the temporal axial gauge 1 one finds 

lUjiu, k) = (s^ - ^) n T (cu, k) - ^n L (cu, k), 
n 00 (u, k) = -u L (cu, k), n 0i (uk) = -^n L (u, k), (2.2) 

where the longitudinal and transversal parts of the self-energy are given by 



IIx, (a;, k) = m 2 D 
and 



u ui + k 
— In - 

2k to — k 



n T (uj,k) = - 



m D -} 



k 2 



-n. 



m 



D 



Ana* 



2N + Nf T2+ N f ^ 



6 



2tt 2 



, (2.3) 



(2.4) 



is the QCD Debye mass squared. In the temporal axial gauge the gluon propa- 
gator then reads 



where 



G L (u,k) = 



kikj 



-1 



G T {u,k) + ^^=G L {u,k), 
k z ui A 



k 2 + U L (u;,ky 



G T (uj,k) = 



-1 



to 2 - k 2 -U T (u,k) 



(2.5) 



(2.6) 



The poles u of the gluon propagator then fall into two branches, correspond- 
ing to a longitudinal (Gj^(ui,{k), k) = 0) and a transversal (G^ 1 (cuxik) , k) = 0) 
mode. The dispersion relations E L (k), E T (k) and associated damping rates 
7l(&),7t(&) are simply given by the real and imaginary part of the poles, re- 
spectively, 

uo L {k) = E L {k) - i lL (k), u T {k) = E T {k) - i lT (k), (2.7) 

or more specific for the above modes, 

= k 2 + ReU T (E L (k), k), E^{k) = k 2 + ReU T (E T (k), k), 
= lmU L (E L {k) + i lL (k), k), lT {k) = -^Imn T (£ T (i) + i lT (k), fc).(2.8) 

In Fig. 2.1, a plot of the dispersion relations in quadratic scale is shown: there 
are propagating modes above a common plasma frequency oo > ui v \ = a/ 
which for the transverse branch tend to a mass hyperboloid with asymptotic mass 
squared 



™2 m D 



(2.9) 



x This particular gauge is chosen for convenience only; a proof of the gauge independence of 
the dispersion relations below can be found in [66]. 



2.1. THE GLUON PROPAGATOR IN THE HTL APPROXIMATION 13 



1.75 


v — * 1 


1 . 5 




1.25 




LU 2 /m 2 D i 


^ ' 

^ ' 


0.75 




. 5 


* 

* 

* 

* 



Figure 2.1: Dispersion relations for the transversal (Gt) and longitudinal (G\l) 
part of the gluon propagator in quadratic scales. Positive values of k 2 correspond 
to propagating modes (mass hyperboloids show as straight lines parallel to the 
light cone indicated by the dashed line) whereas negative values show screening 
phenomena. 



whereas the longitudinal branch approaches the light-cone exponentially with 
a residue that is exponentially suppressed. Indeed, this last mode (which is 
sometimes called the longitudinal plasmon) does not have a vacuum analogue 
and disappears from the spectrum at k — > oo. The transverse mode, however, 
represents an in-medium version of the physical gauge boson polarization. 

For uj < ujpi, k turns out to be imaginary which corresponds to the phe- 
nomenon of screening; in QCD, there is magnetic (transversal mode) as well as 
electric (longitudinal mode) screening with a screening length given by l/\k\ as 
long as uj > 0. 

In the static limit uj — 0, the self-energies become 11^(0, k) = m 2 D and 
IFr(0, k) = 0, respectively. Accordingly, whereas the longitudinal branch pre- 
dicts screening in the electrostatic sector with inverse screening length \k\ — mo, 
the solution to the transverse branch, k — 0, signals the absence of magnetostatic 
screening. In fact, in Abelian gauge theories there are strong arguments that the 
absence of a magnetic screening mass is required by gauge invariance of the the- 
ory [67]. However, in non- Abelian gauge theories this absence is not guaranteed 
and indeed lattice QCD simulations do find a screening behavior in the transverse 
sector [68]. Note that also analytically one can find a nontrivial magnetostatic 
behavior if one no longer requires the system to be isotropic; the result for the 
magnetic screening mass, however, is somewhat unexpected and this case will be 
described in more detail in chapter 6. 

To leading order, the imaginary part of the gluon self-energy vanishes for 
uj 2 > k 2 . Accordingly, the damping rates of the quasiparticles are zero to this 
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order of the coupling. However, there is a contribution at order a s , which has 
been calculated by Braaten and PlSARSKl [69] in the zero momentum limit 
after several incomplete attempts by various authors in the 1980's (see [65] for a 
collection of published values) . For nonzero momentum the damping rates have 
an apparent logarithmic singularity that is cut off by the magnetostatic mass [70], 
which has to be calculated non-perturbatively; only the zero momentum result 
seems to be protected from this effect, although there may be different problems 
arising at next-to-leading order [71]. 

Finally, in the region with u 2 < k 2 the imaginary part of the gluon self-energy 
is large (~ a s T 2 ); this corresponds to the possibility of Landau damping, which 
is the dissipation of energy in the plasma [6] (see also chapter 7 where the energy 
loss of a heavy parton is calculated). 



2.2 The fermion propagator in the HTL approxi- 
mation 

The fermion propagator S (suppressing Lorentz indices) fulfills a Dyson-Schwinger 
equation similar to the gluon propagator 

S- 1 = + S, (2.10) 

where E is the fermionic self-energy. The most general form of S compatible with 
the rotational and chiral symmetries is 

E(u, k) = a(u, k)T° + b(u, k)k ■ T, (2.11) 

where T are the Dirac matrices. This can be rewritten as 

r S(cj, k) = S + (cj, fc)A+(k) - S_(w, fc)A_(k), (2.12) 

where S±(c<j, k) = b(u, k) ± a(u>, k) and the spin matrices 

A ± (k) = i^k, A+ + A_ = l, A| = A ± , (2.13) 
A + A_ = A_A + = 0, TrA ± = 2, (2.14) 

project onto spinors whose chirality is equal (A + ), or opposite (A_), to their 
helicity. Using Sq 1 = —ft + S one can decompose the propagator 

S(uj, k)T = A + {u, k)A + + A_(w, fc)A_, (2.15) 

with Ajjj 1 = — [uj =p (k + S ± )]. Accordingly, there are also two fermionic pole 
branches corresponding to A^ 1 = and Al 1 = 0, where the latter one is occasion- 
ally referred to as plasmino branch. Using the result for the fermion self-energy 
in the HTL approximation, 

v. / ,x M 2 ( cu±k, cu + k\ 

^k) = - !-_!„_ , (2.16) 
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Figure 2.2: Dispersion relations for the fermionic quasiparticles in quadratic scale. 
Shown are the regular (+) and plasmino (-) branch, together with the light cone 
(dashed-line). 



where M 2 is the plasma frequency for fermions, 

Ana s (N 2 - 1) / /x 2 
167V V Ti 



^= r:», j [t 2 +^), (2.i7) 



one can proceed to plot the dispersion relations E + {uj), E_{oj) for the fermionic 
modes, shown in Fig. 2.2. 

The plasmino branch shows an interesting minimum of u> at uj/M ~ 0.93, 
k/M ~ 0.41 and approaches the light-cone for large momenta with exponen- 
tially vanishing residue; the regular branch tends to a mass hyperboloid with 
asymptotic mass squared 

= 2M 2 . (2.18) 

The damping rates also vanish to leading order but there is a contribution to 
order a s , as was the case for the gluonic modes [70]. 



2.3 Simple quasiparticle model for the EOS 

The main building block one uses in the derivation of the HTL approximation 
is that the principal contribution for the loop integral in e.g. the self-energies 
comes from momenta which are of order k ~ T, dubbed "hard" (hence the name 
hard-thermal-loop). Therefore, for weak coupling a s when momenta of order 

1 /2 

k ~ a s T (dubbed "soft") are much smaller than hard momenta, the thermody- 
namic behavior of the system is dominated by the hard excitations. However, 
the relevant collective modes at large momenta have already been identified in 
the previous sections: while the longitudinal plasmon as well as the plasmino 
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mode are exponentially suppressed, the transverse gluonic and positive fermionic 
branches represent particle excitations that propagate predominantly on simple 
mass shells with dispersion relations 

E T (k) 2 = ml + k 2 , E + (k) 2 = M 2 OQ + k 2 , (2.19) 

with the asymptotic masses and given above. It is therefore straightfor- 
ward to model the thermodynamic pressure as a gas of weakly interacting massive 
particles with residual mean-field interaction B [72, 44], 

p(T, fi) = p g (T, ml) + p q (T, fi, Ml) - fl(rr&, M 2 J, (2.20) 

where 

d 3 k 



Pg (T,ml) = -2(^-l)|A^ln(l-e-^), 



(2.21) 



/d 6 k 
^ In (1 + e-^~^ T ) + (jm - -r) (2.22) 

represent the partial pressure from gluons (g) and quarks (q). Using the sta- 
tionarity of the thermodynamic potential under variation of the self-energies one 
finds 

dB dp g (T,ml) dB dp q (T^,M 2 J_^ 



dml dml ' dMl dMl 



implying 



Q _dp_ d[p g (T,ml)+ Pq (T,^Ml)} n _dp_ d Pq (T,fi,Ml) 
S ~dT~ dT ' n ~dfi~ dfi 

for the entropy and number density, s and n, respectively. 

If one expands Eq.(2.20) in powers of the coupling a s , the leading order pertur- 
bative result P0+P2 given in Eqs. (1.2,1.3), and part of the higher order correction 
is reproduced. Unexpanded, however, Eq.(2.20) represents a thermodynamically 
consistent resummation of terms in all orders of a s , suggesting also a possible 
applicability of the model in the strong coupling regime. Indeed, by introducing 
an ansatz for the strong coupling (inspired by the running coupling Eq.(1.4)) that 
contains two phenomenological parameters A and T s at fj, — 0, 



12tt 

[IIN -2N f y^\\(T + T s )/T c X 



(Xs,eS(T,/J, = 0) = 7777; r-\ t\ r \ 1 [ \ / m 1 m \ / m "1 2 ^ (2.25) 



it turns out that one can accurately describe lattice data for the entropy and 
eventually the pressure (see chapter 3 for details). 

Moreover, for non-vanishing chemical potential, one obtains a flow equation 
for the coupling: since the quasiparticle masses depend on T and \i explicitly 
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as well as implicitly through the coupling a s (T,fi), Maxwell's relation ds/dfi = 
dn/dT implies a partial differential equation for a s which takes the form 

aT w + a ^ = b > (2 - 26) 

where the coefficients ar, and b depend on T, \i and a s . Given a valid boundary 
condition, a solution for a s (T, jj) can be found by solving the above flow equation 
by the method of characteristics; once a s is thus known in the T, \i plane, the 
quasiparticle partial pressures p g and p q are also fixed completely. Finally, the 
residual interaction B is then given by the integral 

l=g,q 

where B is an integration constant that has to be fixed either through lattice 
data or through matching the pressure on to a low-energy effective theory. In 
chapter 3 the details of this calculation and also an interpretation of the results 
will be given. 

2.3.1 Properties of the flow equation 

When calculating the coefficients ax, and b of Eq.(2.26), one finds that in 
particular ar and obey [45] 

a T {T, n -> 0) = 0, a M (T -> 0, fi) = 0, (2.28) 

while the respective other coefficient does not vanish in these limits. The char- 
acteristic equations following from Eq.(2.26), 

a^dT = ardfj,, a^da s = bd/i, (2.29) 

then imply that the characteristics are perpendicular to both the T and the \x 
axes. Therefore, specifying the coupling a s on some interval on the T or /x axes 
represents a valid boundary condition for the flow equation (2.26). 

In the limit of vanishing coupling a s — > it can be shown that the coefficient 
b vanishes. The coupling is then constant along the characteristics, which become 
ellipses in the variables T 2 and /i 2 [44], given by 

AN + hN f . , „ 



9N f 



n z \7r/ 



Note that for nonzero coupling this is still true approximately. 

However, the solution to the flow equation is not forcibly unique everywhere, 
since the characteristics may turn out to be intersecting for certain regions in the 
T, /i-plane. This ambiguity as well as its resolution will be treated in more detail 
in chapter 3. 
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2.4 Summary 

In this chapter I have defined and characterized the leading-order gluonic and 
fermionic quasiparticle excitations relevant for QCD. From the respective self- 
energies the dispersion relations were calculated and it was found that for hard 
momenta two of these excitations vanish exponentially while the remaining two 
propagate approximately on simple mass shells. Following Gorenstein and 
Peshier [72, 44], I argued that the thermodynamic pressure of QCD can be 
modeled by that of a free gas of massive particles together with some residual 
interaction, which itself is determined by requiring stationarity of the pressure 
with respect to the self-energies. Anticipating that by using a simple 2-parameter 
ansatz for an effective strong coupling this model can be used to accurately de- 
scribe lattice data for the entropy and pressure (which will be done in chapter 
3), I showed that by using Maxwell's relation one obtains a flow equation for the 
coupling, which, once solved, will give an equation of state of deconfined QCD for 
arbitrary temperature and chemical potential. Finally, I derived the equations 
for the characteristics of the coupling flow equation and anticipated that they 
will resemble ellipses in the variables T 2 and /x 2 , with endpoints perpendicular to 
the T and \i axes. 



Chapter 3 

The HTL quasiparticle model 



In the last chapter I have set up a simple quasiparticle model of the QCD pressure 
that only takes into account the hard excitations, ignoring longitudinal plasmon 
and plasmino quasiparticles. A comparison between the strictly perturbative 

5/2 

pressure known to order a s and a perturbative expansion of the simple quasi- 
particle model shows that - while the leading orders are correctly reproduced 
- only l/A/y/2 or about 18% of the order a^J 2 term are included in the sim- 
ple quasiparticle model. However, it is precisely the inclusion of this term (also 
called the plasmon term) that makes the result for the pressure go "wild" in the 
first place, spoiling the convergence of the strictly perturbative expansion as was 
shown in Fig. 1.2. Consequently, neither does the simple quasiparticle model 
incorporate a big part of the plasmon effect (which could mean quantitative re- 
sults are incorrect by an unknown amount) nor is it a priori clear whether an 
extension including the full plasmon term would not go astray as violently as 
in the strict perturbative case (implying that also qualitatively the results from 
the simple quasiparticle model would be dubious). Therefore, it is important to 
consider refinements of the simple quasiparticle model that incorporate more of 
the physically important plasmon term, so that the above caveats can be settled. 

One such refinement, namely that of a quasiparticle model for the QCD ther- 
modynamic pressure that is based on the HTL-resummed entropy [28, 73, 29], 
will be considered in the following. 

3.1 The entropy of the QCD plasma 

As a starting point, one uses the result that the thermodynamic potential Q = 
—pV can be written as the following functional of the full gluon and fermion 
propagators G and S [28]: 

-tt[G,S} = ^TrlnCr 1 -Trln^ 1 - -TyIIG + TrZS + <5>[G, S], (3.1) 
_£ 2 2 
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Figure 3.1: Diagrams for <3> at 2-loop order in QCD in the temporal axial gauge. 
Wiggly lines refer to gluon, plain lines refer to fermions. 



where IT and E are the gluon and quark self-energies of chapter 2 (related to the 
full propagators by Dyson's equation (2.1, 2.10)) and Tr includes traces over color 
indices and also over Lorentz and spinor indices when applicable. &[G, S] is the 
sum of the 2-particle-irreducible "skeleton" diagrams, which to 2-loop order are 
shown in Fig. 3.1. 

The essential property of the functional Q[G, S] is to be stationary under 
variations of G and S (at fixed G and So), 

jg=0, ^ = 0, (3.2) 

so the physical pressure is then obtained as the value of fl at its extremum. Note 
that the stationarity conditions imply 

which together with the Dyson's equations define the physical propagators and 
self-energies in a self-consistent way. Accordingly, by selecting a class of skeletons 
in &[G, S] and calculating II and S from the above relations one obtains a self- 
consistent approximation of the physical pressure. 

Using the temporal axial gauge which preserves rotational invariance and 
also makes the Faddeev-Popov ghosts decouple from the theory, the gluon and 
quark propagators can be decomposed as in the previous chapter, so that one can 
evaluate Eq.(3.1) using standard contour integration techniques. Accordingly, the 
entropy density s = — dn J^ is evaluated to be 



Sg + s g + s', (3.4) 



where 



s 9 = -(iV 2 -I) J {2 [Imln (-c 2 + k 2 + U T ) - lmU T ReD T ] 

+Imln (k 2 + n L ) + Imn L ReD L } , 
s, = -2AW / | I g I ^{ImlnA;' + Imln(-A:') 
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s = — 



-Im£ + ReA+ + Im£_ReA+} , 
d(T/V$[G,S}) r ^ l 



V$[G,S\) /Ar2 . f d A k dn(u) . n „ T ^ „ „ T ^ , 
dT [ U d + (N 2 -l) J j^i-^j 1 (2Ren T ImD T -Ren L ImD i ) 

+2iViV/ / (2k)* } { ReS + ImA + - ReS_ImA + } . (3.5) 

Here ra(u;) and f(uj) are the Bose-Einstein and Fermi-Dirac distributions, 

1 1 1 

U(U}) = expu/T-V /M = exp ( w - /x)/T + 1 + exp (u + fj)/T + 1' (3 ' 6) 

and decomposition of the propagators from chapter 2 has been used. 

When restricting to an approximation where $[-D, S] contains only 2-loop 
skeletons (shown in Fig. 3.1), one finds the remarkable property that 

s' = 0, (3.7) 

which has been shown in Refs. [74, 29] for QED and QCD, respectively. Since 
in this approximation the self-energies and propagators are to be determined 
self-consistently by solving "gap-equations" of the form 

D T X = -u 2 + k 2 + Tl T [D T , D L , A+, A_], (3.8) 

where Ilr is given by one-loop self-energy expressions only, the resulting expres- 
sion for the entropy (3.4) is also effectively one-loop. 

This has important consequences: firstly, one has reduced the problem of 
calculating a two-loop quantity (the free energy or thermodynamic pressure) to 
a simple integration of an effectively one-loop quantity (the entropy expression 
from above); secondly, since the term s' has been interpreted as the contribution 
from the residual interaction of the quasiparticles [29], the result s f — means 
that this contribution vanishes for the QCD entropy in the 2-loop approximation. 
I will make use of this interpretation in the following. 



3.1.1 Approximately self-consistent solutions 

Unfortunately, "gap equations" like Eq.(3.8) are non-local, which makes their 
exact solution prohibitively difficult. Moreover, the gap equations contain UV 
divergencies and therefore require a renormalization scheme compatible with the 
self-consistent structure of this non-local equation. Whereas this daunting task 
has been at least formally resolved in the scalar A0 4 model [33, 34], the richer 
structure of QCD along with unresolved problems concerning gauge-dependencies 
of the two-particle-irreducible scheme at higher orders [75] render a self-consistent 
renormalization of the gap-equations in QCD an open problem. 

Therefore, it is convenient to construct so-called "approximately self-consistent" 
solutions, which maintain the properties of Eq.(3.3) together with the Dyson's 
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equations up to and including a given order of a s , and which are manifestly 
gauge independent and UV finite. For example, using the HTL results for the 
self-energies and the resulting full propagators (given in chapter 2) constitutes 
just one such approximately self-consistent solution; therefore, using the HTL 
approximation for the self-energies and propagators in Eqs. (3.4,3.5,3.7) gives a 
gauge-invariant, non-perturbative approximation to the full entropy. 



3.2 The HTL quasiparticle model 

Similar to the simple quasiparticle model, one can construct a model for the 
thermodynamic pressure that incorporates the relevant HTL quasiparticle exci- 
tations, which I will call the HTL quasiparticle model in the following [76]. Since 
this model should be a refinement of the simple quasiparticle model, two require- 
ments have to be fulfilled: firstly, the HTL model has to incorporate more of the 
physically important plasmon effect, and secondly, it should constitute a natural 
extension of the simple quasiparticle model, including its structure in some sense. 

The key to the construction of the HTL quasiparticle model lies in the obser- 
vation that similar to the results in the simple quasiparticle model the entropy in 
the two-loop ^-derivable approximation of the previous section is given in terms 
of HTL quasiparticle excitations only, while the contribution from residual inter- 
actions vanishes. One therefore defines the HTL quasiparticle model entropy to 
be given by Eq.(3.5), which is known to be a fairly good approximation to the 
full entropy of the quark-gluon plasma [29]. Taking all this into account, one is 
led to the following ansatz for the thermodynamic QCD pressure: 

p(T, fi) = Pg ,HT L (T, m 2 D ) + P(1 ,htl(T, fi, M 2 ) - B HTL (m 2 D , M 2 ), (3.9) 

with 

P 9 ,htl = ~2(N 2 " 1) / ^) [2Imln (-c 2 + k 2 + H T ) 

-2Imn T ReD T + Im In (k 2 + n L ) + Imn L ReD L ] (3.10) 
— j( -/H[Imln(k-u; + E + ) 

-Im£ + ReA + + Im In (k + u + E_) + Im£„ReA_] , (3.11) 

obeying the stationarity conditions 

dp dp 



dm 2 D dM 2 



0. (3.12) 



Clearly, the model (3.9) fulfills all of the above requirements: the pressure is de- 
scribed by a gas of weakly interacting HTL quasiparticles, which have momentum- 
dependent dispersion relations that were given in chapter 2. Therefore, the in- 
tegrations in Eq.(3.11) pick up contributions equal to those of the simple quasi- 
particle model for large momenta, while for soft momenta the richer structure of 
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Figure 3.2: Entropy data generated from Ref. [18] vs. fitted model entropy. The 
HTL and simple quasiparticle model fits are nearly indistinguishable to the naked 
eye. 



the quasiparticle excitations are taken into account. Consequently, when doing 
a perturbative expansion of Eq.(3.11) at fi = in powers of the coupling a s , one 
finds that 25% of the plasmon effect is included in the HTL model, in contrast 
to 18% for the simple quasiparticle model. Furthermore, the stationarity condi- 
tion ensures that there is no residual interaction to the entropy (s' = 0) and the 
resulting model entropy is indeed given by s = s g + Sf. 

In its setup, the HTL quasiparticle model is very similar to the simple quasi- 
particle model, so the general discussion of this model in chapter 2 remains valid 
to a large extent also for the HTL model. Indeed, to tell whether the plasmon 
effect has a destabilizing effect in these quasiparticle models similar to what hap- 
pens in strict perturbation theory, one has to compare the results from both 
models in a standardized procedure, so I will follow the steps from Eq.(2.25) 
onwards for both models, filling in the details and pointing out the differences. 

3.3 Modeling the lattice entropy and pressure 

Restricting /i = and using the ansatz (2.25) for the strong coupling a s , the 
entropy of the HTL and simple quasiparticle models becomes a function of the 
temperature T and two fit parameters T s and A, s = s(T,T s ,\). Using lattice 
data for Nf = 2 and N = 3 from Ref. [18] together with an estimated continuum 
extrapolation from [45], these two parameters T s and A are determined by a least 
square fit to the data. The numerical values for these fit parameters (which are 
reproduced in chapter 4) turn out to be quite similar, and the respective fits to 
the lattice data (shown in Fig. 3.2) are nearly indistinguishable. 
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Figure 3.3: Residual interaction B for the HTL (full line) and simple model 
(dashed-line), respectively. 

Once the fit parameters have been fixed, so has a s ^ at vanishing chemical 
potential, and consequently also the quasiparticle masses. One can then use 
Eq.(2.27) together with dfi = 0, 

*- (3 - 13) 

i=g,q 

to determine the residual interaction as a function of the temperature up to an 
integration constant B . This integration constant has been determined by a fit 
of the lattice pressure with the simple quasiparticle model in [45], obtaining 

B = B\ Tc = 1.1 T 4 , (3.14) 

so that the resulting model pressure at T c becomes 

p(T e ,n = 0) = 0.536(1)T C 4 . (3.15) 

To facilitate comparison, I fix the HTL model value of B by requiring 
Phtl(T c , H = 0) ~ 0.536 T 4 , obtaining B HTL = 0.82 T c 4 . The resulting residual 
interactions for the HTL and simple quasiparticle models are plotted in Fig. 3.3; 
finally, the overall model pressure is compared to the lattice data in Fig. 3.4. 

Clearly, both the HTL and the simple quasiparticle model can be used to 
describe the lattice data very accurately once a 2-parameter fit for the effective 
coupling is adopted, even close to the deconfinement transition. Furthermore, 
since the fit parameters turn out to be only slightly different and the resulting 
fits to the lattice data are nearly indistinguishable, one concludes that (at least 
at p, — 0) the inclusion of the plasmon effect does not seem to destabilize the 
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quasiparticle models; on the contrary, even though the effective coupling gets 
rather large near T = T c (see chapter 4 for details), the effects of changing the 
percentage of the included plasmon effect are rather small. Therefore, one gains 
some confidence that a phenomenological description of the thermodynamic QCD 
pressure through quasiparticle models is also applicable in the non-perturbative 
regime, and in turn may be used to derive results that lie parametrically far from 
the initial (lattice) data. 



3.4 Solving the flow equation 

At non-vanishing chemical potential, the models are required to fulfill Maxwell's 
relation 

ds dn . . 

m = w <3 - 16) 

Since the entropy s as well as the number density n depend on T and \i both 
explicitly through the Bose-Einstein and Fermi-Dirac distribution functions as 
well as implicitly through the quasiparticle masses, one has 

ds dn sr-^ ds dm 2 dn dm 2 

111 ~ dr = 2s drf~d^ ~ dn^lr' ^ 3 ' 17 ' 

i=g,q 

because the explicit derivations of s and n with respect to // and T cancel, re- 
spectively. The quasiparticle masses themselves also depend on T and /i both 
explicitly and implicitly through the coupling a s , e &, so that e.g. for the Debye 
mass one finds 
dm 2 . 



dfi 



f = 5#0*™.<T,rt + ^^4, (^T> + (3.18) 
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Figure 3.5: Comparison of the shape of the characteristics for the HTL (full lines) 
and simple quasiparticle model (dashed lines). T c denotes T c |/x=o- 



Therefore, Eq.(3.17) represents a partial differential equation for the coupling 
that can be written in the form already anticipated in Eq.(2.26), 



da s da s 



b; 



(3.19) 



since it takes a considerable amount of algebra to derive the coefficients Or, 
and b in the HTL model, their explicit form may be of use in further studies. 
However, since their appearance is rather unwieldy, I refrain from reproducing 
these coefficients in the main text and have relegated them to appendix B. 

When evaluating the coefficients numerically and comparing their value for 
the HTL and simple quasiparticle model, one finds that while and a M assume 
nearly equal values, the coefficient b turns out to differ noticeably which has 
interesting consequences, as will become clear in the following. The solution to the 
flow equation (2.26) for the coupling is found by using the initial value condition 
a s,es(T,fi = 0) and then numerically solving the characteristics of Eq.(2.26), 



a^dT 



Qj ^ d Gt g 



bdfi. 



(3.20) 



The shapes of the characteristics are shown in Fig. 3.5 for both the HTL and the 
simple quasiparticle model; as anticipated in chapter 2, the characteristic curves 
resemble ellipses in T 2 and /z 2 , which are perpendicular to both the T and the 
jj axis. Furthermore, for the simple quasiparticle model there is a set of char- 
acteristics (those which start at \i = in the interval [T c , 1.06TJ) that intersect 
in a narrow region, indicating that there the solution to the flow equation is not 
unique (see Fig. 3.6). However, when also calculating the pressure along the 
characteristics it turns out that it becomes negative in the intersection region; 
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Figure 3.6: Characteristics starting very near T c : for the simple quasiparticle 
model (dashed-lines) intersections are occurring, while for the HTL model (full 
lines) this is not the case. 



it has been argued [45] that this implies that a transition to another phase with 
positive pressure has been occurring somewhere outside this region, which there- 
fore means that the quasiparticle description had broken down and the ambiguity 
would be of no physical relevance. 

Interestingly, for the HTL model, the difference in the flow equation param- 
eter b seems to completely lift the ambiguity problem: the characteristic curves 
are well behaved and do not intersect anywhere, as can be seen in Fig. 3.6. Un- 
fortunately, the pressure still turns negative for some region in the T, /i plane 
for the HTL model, so that clearly also this model breaks down for certain T, \i. 
On the other hand, one can turn this breakdown of the quasiparticle models 
into a virtue by assuming that the deconfinement transition occurs near the line 
T = T(/i) where the pressure vanishes. Therefore, an estimate for the "critical" 
transition line for arbitrary \i becomes calculable within the quasiparticle models, 
(see chapter 4 for details). 

3.5 The pressure at finite chemical potential 

Solving the flow equation provides one with an effective coupling in the whole 
T, n plane, so that - via the explicit form of the quasiparticle masses - the 
contributions p g and p q to the full pressure (2.20,3.9) are fixed. It remains to 
calculate the residual interaction £>, which is most easily done by integrating 
Eq.(2.27) along the characteristic curves of the flow equation (2.26). The result 
for the full pressure is shown in a 3D-plot in Fig. 3.7. 

As will be discussed in more detail in chapter 4, this result is in fairly good 
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Figure 3.7: The full pressure of the HTL quasiparticle model normalized to the 
free pressure po in the T, \i plane. 
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agreement with various lattice approaches. On the other hand, Ipp and Rebhan 
[77] have found that in the exactly solvable large- Nf model of QCD the pressure 
has a conspicuous "kink" near \i ~ nT, signalling an abrupt change of the pressure 
for larger chemical potentials. This kind of behavior is not reproduced in the 
quasiparticle model approach for 2 flavors; on the contrary, it turns out that 
the quasiparticle model pressure is well described by its small-/i behavior even 
for large values of the chemical potential, as will be discussed in more detail in 
chapter 4. However, it is by no means clear that results obtained in the large- Nf 
case are also relevant in general for the physically interesting situations of Nf = 2 
or Nf = 3, although for special cases their effect may well be considerable, as has 
recently been demonstrated [78]. 

3.6 Summary 

In this chapter I have introduced the entropy of the quark gluon plasma through 
the so-called ^-derivable formalism. A 2-loop approximation for the functional 
$ that allows the entropy to be written as an effectively one-loop quantity was 
presented, and it was pointed out that because of problems concerning renormal- 
ization and gauge-dependencies it is useful to introduce the notion of approximate 
self- consistency. Since using the HTL self-energies together with the respective 
Dyson's equations for the propagators introduced in chapter 2 corresponds to such 
an approximately self-consistent approximation to the full quark-gluon plasma en- 
tropy, this expression was used as a basis for an HTL quasiparticle model, and it 
was shown that this model corresponds to an extension of the simple quasiparticle 
model introduced in chapter 2. 

Fitting two flavor lattice data for the pressure it was found that at vanish- 
ing chemical potential the HTL model constitutes only a small correction to the 
simple quasiparticle model, so it seems that the plasmon effect which plays such 
a devastating role in strict perturbation theory is not destabilizing the quasi- 
particle implementation. Moreover, an explicit solution to the flow equation for 
the coupling shows that while the simple quasiparticle model solution becomes 
ambiguous in some region in the T, \x plane, one does not encounter any such 
problems for the HTL model. 

Finally, I investigated the pressure at finite chemical potential and discussed 
its general behavior in view of recent exact results for the large-flavor limit of 
QCD. 
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Chapter 4 



The QCD EOS for arbitrary fi and 
completion of plasmon effect 

Having demonstrated in the last chapter that lattice data at vanishing chemical 
potential can be accurately described by quasiparticle models with an effective 
coupling, I now concentrate on the quantitative predictions one can make by 
mapping lattice data to nonzero chemical potential. More precisely, first the 
quasiparticle model results at low chemical potential will be compared to those 
obtained through lattice calculations and strictly perturbative methods, and then 
the case of high chemical potential and low temperatures, where nearly all other 
methods eventually break down, will be investigated. However, I will ignore the 
effect of color superconductivity, which should become important at very small 
temperatures only, assuming that because of the comparatively small energy gap 
it constitutes only a minor modification to the thermodynamic pressure in the 
region of interest [46]. The results in this chapter are based on Ref. [79]; however, 
here I will be able to go into more detail and also add some new results. 

In the previous chapters, I have introduced the HTL quasiparticle model that 
includes 1/4 of the full plasmon effect at vanishing chemical potential, which is 
a factor of a/2 more than in the simple quasiparticle model that uses only the 
asymptotic HTL masses. Motivated by the rather successful generalization of the 
simple to the HTL quasiparticle model I will try to improve the accuracy of the 
HTL quasiparticle model further by setting up a model that includes even more 
of the perturbative plasmon effect. In principle, for this task one would need the 
complete (momentum-dependent) next-to-leading order corrections to the self- 
energies, which up to now no one had the strength to calculate; however, it has 
been shown in Ref. [28, 29] that - as far as the plasmon effect is concerned - the 
contribution to the thermodynamic quantities of these next-to-leading order cor- 
rections may be approximated by their averaged contribution to the asymptotic 
masses. Therefore, this opens the possibility of setting up a model that con- 
tains the full plasmon effect at /i — 0, albeit at the price of introducing another 
unknown parameter, as will be discussed in the following. 
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4.1 NLO quasiparticle models 

Denoting the momentum-dependent next-to-leading {o?J 2 ) order corrections to 
the self-energies as 5Ut(oo, k) and 5T> + (the other branches do not receive contri- 
butions at this order), respectively, the averaged corrections to the asymptotic 
masses (2.9,2.18) have been calculated to be [28, 29] 

- 9 f dk kn'(k)Re5HT(uj = k) 

5m ~ = 1 — jskw — 1 = ~ 2a - NTm ° (41) 

and similarly 

- 2 [dkk(fUk) + f(k))Re2k5X + (u = k) (N 2 -l) , x 

For the values of the coupling a s considered here, these corrections are so large 
that they give tachyonic masses when treated strictly perturbatively. In Ref. [29] 
it has been proposed to incorporate these corrections through a quadratic gap 
equation which works well as an approximation in the exactly solvable scalar 
0(N — > oo)-model, where strict perturbation theory would lead to identical dif- 
ficulties. However, for the fermionic asymptotic masses, in order to have the 
correct scaling of Casimir factors in the exactly solvable large- A^j limit of QCD 
[80], a corresponding gap equation has to remain linear in the fermionic mass 
squared. Choosing the correction therein to be determined by the solution to the 
gluonic gap leads to [81, 30] 

™L = - 2V2a s NTm oc (4.3) 
»-> 2 ,, 2 V2a s (N 2 - 1)T 

Ml = M{ y -'"x- (4.4) 

where m 2 ^ and are the leading-order gluonic and fermionic asymptotic masses 
as given in (2.9) and (2.18). This in fact avoids tachyonic masses for the fermions 
as long as Nf < 3. 1 

Finally, since the averaged quantities fh 2 ^ and are the effective masses at 
hard momenta only, a cutoff scale A = ^/2iiTm D c^ is introduced that separates 
soft from hard momenta. The quasiparticle pressure for this model, which in the 
following will be referred to as NLA-model, then separates into a soft and a hard 
component for both gluons and fermions. The soft contributions are given by 
expressions similar to Eq.(3.11), but with A as upper limit for the momentum 



1 For Nf = 3 the solutions to the above approximate gap equations happen to coincide with 
those obtained in the original version of two independent quadratic gap equations of Ref. [29]. 
For the case Nf — 2 considered here, the differences are fairly small. For Nf > 3, however, the 
necessity to avoid tachyonic masses would restrict the range of permissible coupling strength 
a*. 
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integration. For the hard contributions, the momentum integrations run from A 
to oo and the mass pre-factors in the HTL self-energies (2.3,2.16) are replaced by 
their asymptotic counterparts, m 2 D — > Zm 2 ^ and M 2 — > \M^ Q . 

The single free parameter c\ in A can be varied around 1 to obtain an idea of 
the "theoretical error" of the model. In the following the range ca = \ to ca = 4 
will be considered; note that ca = oo corresponds to the HTL-model (which has 
been used as cross-check) since all hard corrections are ignored. On the other 
hand, c\ = would assume that (4.3) and (4.4) represent good approximations 
for the NLO corrections to the spectral properties of soft excitations. However, 
the few existing results, in particular on NLO corrections to the Debye mass [82] 
and the plasma frequency [83], appear to be rather different so that it seems safer 
to leave the soft sector unchanged by keeping a finite c\. 

4.2 Fitting lattice data at fi = 

Going through the program of fitting the entropy expressions from the models 
under consideration to lattice data [18] for Nf = 2 as described in chapter 3, one 
finds the following values for the fit parameters T s and A: 





simple 


HTL 


c A = 4 


c A = 1 


c A = 1/4 


T s /T c 


-0.89 


-0.89 


-0.89 


-0.84 


-0.61 


A 


17.1 


19.4 


18.64 


11.43 


3.43 



It can be seen that the results for the simple, HTL and NLA quasiparticle model 
with ca = 4 are very close. The fits to the entropy data all lie in a narrow band 
for all the models considered [79], resembling the fits for the simple and HTL 
quasiparticle models shown in Fig. 3.2. The fitted effective coupling a s for the 
various models is shown in Fig. 4.1; for comparison, also the 2-loop perturbative 
running coupling in MS is shown, where the renormalization scale is varied be- 
tween 7rT and 4nT and T c = 0.49Aj^g. As can be seen from the plot, the results 
for the effective coupling are well within the range of the 2-loop perturbative 
running coupling (for the case ca = 1/4 and renormalization scale ttT the results 
even seem to be identical, which is, however, probably only a coincidence). The 
result for the coupling obtained in the semi-classical approach of Ref. [84] 2 is also 
shown in Fig. 4.1. 

In general, one can see that the effective coupling becomes bigger when ca 
gets smaller; this is because the hard masses (for equal values of the coupling) are 
smaller than the soft masses, which makes the entropy increase when the hard 
parts become more important. Accordingly, the coupling has to rise in order for 

2 For finite /i and constant temperatures, however, the coupling obtained in [84] rises while 
the quasiparticle model results indicate a decrease of the coupling (which is consistent with the 
standard QCD running coupling with renormalization scale proportional to \JT 2 + (fi/n) 2 ). 
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T/T c 

Figure 4.1: Effective coupling: NLA model results for ca = 4,1 and 1/4 (full lines 
from bottom to top), 2-loop perturbative coupling in MS (gray band) and result 
from Ref. [84] (dashed line). 



the entropy to match the data (therefore, for the extreme case ca = one finds 
huge values of the effective coupling constant). 

The fact that the fitted effective coupling agrees so well with the band of the 
2-loop perturbative running coupling suggests that also simply using the pertur- 
bative running coupling as an input might lead to results for the thermodynamic 
pressure that are not very different; this approach will be tested in chapter 5. 

Once the effective coupling is known for fi = one can proceed to fix the 
remaining integration constant B as indicated in chapter 3. Setting p(T c ) = 
0.536(1) T c 4 one finds 





simple 


HTL 


c A = 4 


c A = 1 


c A = 1/4 


#o/T c 4 


1.1 


0.82 


0.73 


0.6 


0.47 



The fits for the model pressure for the various models also turn out to lie in 
a narrow band [79] resembling the HTL and simple quasiparticle model result 
shown in Fig. 3.4. Therefore, also the NLA quasiparticle model can be used to 
accurately describe lattice data for the entropy and pressure at vanishing chemical 
potential. In contrast to what was found from the comparison between the HTL 
and simple quasiparticle model, the fit parameters indicate that at least for c\ < 1 
the NLA model differs considerably from the other models. Unfortunately, it is 
hard to say if this difference is due to the full inclusion of the plasmon effect in 
the NLA model and thereby represents a better approximation of the underlying 
physics than the HTL model, or rather to the fact that only the averaged NLO 
corrections to the self-energies were taken into account and the deviations of 
this model signal the onset of the breakdown of applicability, which probably 
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occurs 3 at c\ <C 1. In any case, the results for physical quantities from the NLA 
models with 1/4 < ca < 4 turn out to fulfill the requirement that they represent 
again small corrections to the HTL and simple quasiparticle model results, giving 
further confidence that the stability of quasiparticle models is unaffected by the 
plasmon effect. 



4.3 Results for small chemical potential 

The NLA model flow equations for the effective coupling can be solved along the 
lines of the HTL and simple model in chapter 3; the shape of the characteristics 
resembles those of the other models and the crossing of characteristics found in 
the simple quasiparticle model also does not occur in the NLA models (similar 
to what has been found for the HTL-model). 



4.3.1 Susceptibilities 

Once the flow equations of the various models have been solved, it is straightfor- 
ward to calculate the quark- number susceptibilities at /i — 0, 



X(T) 



(4.5) 



The result - normalized to the tree level result xo = ~i>rT 2 - is compared to 
lattice data for N f = 2 [85] and N f = 2 + 1 [40] in Fig. 4.2. As can be seen, there 
is a very good agreement between the model predictions and the lattice data 
(which is completely independent from the one used in determining the model fit 
parameters). The result from a strictly perturbative calculation [86] (not shown) 
also is consistent with the results of Fig. 4.2 at larger temperatures, although 
the uncertainty in the latter is much bigger than that of the quasiparticle model 
results (see also chapter 5). 

Calculating the pressure at finite chemical potential p(T, //) as described in 
chapter 3, one is able to numerically extract the higher order susceptibility 



(4.6) 



by doing a least-square fit to the model pressure 



p m (T, = p(T, 0) + + (4.7) 



3 This hypothesis is strengthened by the fact that although the NLA model with ca = can 
be made to describe the lattice data at [i = 0, the values of the fit parameters for this model do 
not resemble anymore those of the other quasiparticle models given above, and its subsequent 
extension to finite chemical potential gives results that are at odds with respect to the other 
models as well as lattice calculations. 
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Figure 4.2: Susceptibilities from NLA 
models (light gray band), from lattice 
data for Nf = 2 [85] (triangles) and for 
2+1 flavors [40] (boxes). 



Figure 4.3: The tree-level result 
(dotted line) for the quantity x(T) 
compared to the HTL model result (full 
line), lattice results from Ref. [42] as 
well as results from strict perturbation 
theory [37] (dashed line). 



Parameterizing 



X(T) 4! 
2T 2 x(T) 



(4.8) 



one obtains the result shown in Fig. 4.3. In this figure, also the lattice results 
for the corresponding quantity in Ref. [42] (which however have not been ex- 
trapolated to the continuum limit) and strictly perturbative results [37, 87] (for 
the specific renormalization scale p = 8.112T [88]) are shown, indicating general 
agreement between the different approaches. 

A measure of how well Eq.(4.7) describes the quasiparticle result for the pres- 
sure can be obtained by considering the quantity Sp = p(T,p) — p m (T, p), shown 
in Fig. 4.4: in the region considered (p < 10T C ) the biggest deviations occur 
for temperatures near T c ; at T = T c there is a maximum deviation of 8% near 
p ~ 3.5 T c , while for temperatures T > 1.3 T c the deviations stay below 3 x 1CT 3 . 

In fact, this agrees with the result of Fodor and Katz [40], who found that 
a scaling relation Eq.(4.7) without the p A coefficient works rather well for small 



4.3.2 Lines of constant pressure 

Once the data for the pressure at finite chemical potential has been calculated, 
it is also straightforward to extract lines of constant pressure p(T,p) = p(T ,0). 
Clearly, the line p(T,p) = p(T c ,0) is of special interest, since it is believed that 
close to this line the transition from the hadronic to the quark-gluon plasma phase 
is occurring. For small p, a quantity that can be compared to recent lattice results 
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Figure 4.4: Deviations of the quasiparticle pressure from the model pressure 
Eq.(4.7) as a function of \i scaled with the free pressure p : shown is HTL model 
data for T = T c (triangles), T = 1.1T C (stars) and T — 1.3T C (boxes). For higher 
temperatures, the deviations are even smaller than those for T = 1.3T C . The 
result for T = T c reaches to a maximum of about 8%. 



is the slope of this constant pressure line, which is readily calculated for 

the quasiparticle models: 





simple 


HTL 


c A = 4 


c A = 1 


CA = 1/4 


rp d'l 
±c dfJ? 


-0.0634(1) 


-0.06818(8) 


-0.06810(6) 


-0.06329(34) 


-0.041(9) 



One finds that - as anticipated - the simple and HTL quasiparticle models as 
well as the NLA models for ca > 1 give results that lie rather close to each other 
whereas for ca = 1/4 the NLA model result is already somewhat off the other 
values. In Fig. 4.5 the quasiparticle model results are compared to recent lattice 
data for 2 + 1 flavors [40] and 2 flavors [36]. One can see that the band for the 
2 + 1 flavor lattice calculation (resulting from fitting the data with a second- 
order (upper limit) and fourth-order (lower limit) polynomial in \i as proposed 
in [89]) is in very good agreement with the simple, HTL, and NLA c\ > 1 
calculations, whereas for the NLA c\ = 1/4 model the slope is somewhat flatter. 
The lattice study for 2 flavors (dashed lines in Fig. 4.5) predicts a slope of T c ^- = 
—0.107(22), which is significantly steeper. 

Although there is no perfect match between all the various lattice and quasi- 
particle model results for the constant pressure slope, the deviations are small 
enough that there is overall consistency between the different approaches. Indeed, 
it is plausible that the remaining disagreement is due to the fact that none of the 
lattice data have been rigorously extrapolated to the continuum limit, so that 
both the data and the fits of the quasiparticle models are still likely to change 
somewhat when this will be done eventually. 
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Figure 4.5: Lines of constant pressure: NLA models for Nf = 2 from ca = 4, 
c A = 1 and ca = 1/4 (dotted lines from lowest to highest); lattice data for 
N f = 2 + 1 [40] (light gray band) and Nf = 2 [36] (long dashed-lines). 

4.4 Large chemical potential 

Extending the lines of constant pressure from the different models to very small 
temperatures one obtains a crude estimate of the deconfinement transition line 
in this region of phase space (assuming that color superconductivity has only a 
minor effect, as argued above). Denoting with /i c the chemical potential where 
(for vanishing temperature) the pressure at vanishing temperature equals the 
pressure at n = 0, p(0, fi c ) = p(T c ,0), one finds (assuming T c = 172 MeV for 
easier comparison) 





simple 


HTL 


c A = 4 


ca = 1 


c A = 1/4 




548 MeV 


533 MeV 


536 MeV 


558 MeV 


584 MeV 




525 MeV 


509 MeV 


511 MeV 


537 MeV 


567 MeV 



Here /io denotes the value of /i where the pressure vanishes, which may be taken 
as a definite lower bound for the critical chemical potential within the respective 
models. 

In general, these results are in agreement with the estimates for /z c from 
Refs. [89, 45]; for the QP models considered, the lowest and highest results for /i c 
are obtained for the HTL model and the NLA model with ca = 1/4, respectively, 
while the /i c of the simple QP model lies between the NLA ca = 4 and ca = 1 
model values. However, even the lowest result for \i c turns out to exceed the 
value for the critical chemical potential expected in Ref. [90, 46]. 

The quasiparticle results for the pressure at T/T c ~ 0.01 has been compared 
to the result of the perturbative pressure at vanishing temperature in Fig. 4.6. 
The latter is originally due to Freedman, McLerran and Baluni [52, 51] but 
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Figure 4.6: The pressure at vanish- 
ing T: shown are perturbative results 
(light gray band) and NLA models for 
ca from 1/4 to 4 (dark band). 



Figure 4.7: Effective coupling from 
NLA models with c A = 4 to 1/4 (dark 
band) and perturbative 2-loop running 
coupling (light gray band). 



has been recalculated more precisely by Vuorinen [37], 



p(T = 0,n) 



tSNf 
4tt 2 

+ (H 



.a. 



71 



18- 11^2 - 0.53583^/ + ^ In 



Nta. 



71 



-AMln^ 

3 /; /i 



a 



— +0(a d s \na s 

71 



(4.9) 



where in this [91] form p, again denotes the renormalization scale in the MS 
renormalization scheme. The perturbative result has been evaluated using the 
standard two-loop running coupling Eq.(l.G) with Aj^g = T c /0.49 [13] and renor- 
malization scale varied from /i to 3/x, as has been considered in Ref. [90]. It can 
be seen in Fig. 4.6 that while in general the quasiparticle model and perturbative 
results are consistent, the variation in the latter (due to the uncertainty in the 
renormalization scale ft) is much larger than the narrow band obtained from the 
evaluation of the different quasiparticle models. The perturbative two-loop run- 
ning coupling at vanishing temperature itself is shown in Fig. 4.7 as a function 
of chemical potential, compared to the result of the effective coupling for the 
quasiparticle models. 



4.4.1 EOS for cold deconfined matter 

By calculating the number density at T/T c ~ 0.01 and using the results for the 
pressure one obtains an equation of state for cold deconfined matter. As is the 
case for the simple quasiparticle model [45], one finds that also for the HTL and 
NLA models the energy density e is well fitted by the linear relation 

e(jp) = 41? + ap 



with 
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Figure 4.8: Mass-radius relations of non-rotating quark-stars: shown are the 
results using the EOS from the HTL (red full line) and NLA c\ = 1/4 and c\ = 1 
models (blue dashed lines). Also shown are the inferred mass-radius relations of 
RX J1856. 5-3754 from Ref. [92, 93] (light gray region), a possible candidate for 
a quark star. The dark shaded region is excluded by the condition R > 2GM. 
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The value of a ~ 3.2 for Nj = 2 seems to be model independent, in contrast to 
the bag constant B 1 / 4 , which varies between 314 and 360 MeV. 



4.5 Application 1: Quark stars 

As a first application I will use the above equation of state for cold deconfined 
matter to determine the mass-radius relations of non-rotating quark-stars. If 
these stars exist, they would differ from ordinary neutrons stars by being more 
compact and therefore considerably smaller than current neutron-star models 
allow [46], as I will show in the following. 
Starting from Einstein's equations 

R„u - \g^R = 87iT AB/ (e(p),p) (4.10) 

where R^ u , g^ u and R denote the Ricci and metric tensor and Ricci scalar, respec- 
tively. Using the Schwarzschild solution for the metric of a spherically symmetric 
and static star one readily finds the familiar Tolman-Oppenheimer-Volkoff equa- 
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tions [94, 95, 47] 



dp(r) 



G [e(r) +p(r)\ [M(r) + 47rr 3 p(r)] 



(4.11) 



dr 




r 



where 



M(r) =4:ir e(r)r 2 dr 



(4.12) 



o 



and G is Newton's constant. The equations are easily solved by specifying a 
central pressure or density and calculating the size and mass of the corresponding 
star using the equations of state from the previous section. The results are 
shown in Fig. 4.8, where it can be seen that the quasiparticle model EOS allows 
for quark-stars with masses of ~ O.8M and radii of less than 5 km (which is 
similar to the results found in Refs. [46, 45, 96]). Current neutron star models 
- which are consistent with the pulsar data [46] - predict a maximum mass in 
the range of 1.4 to 2 solar masses with radii ranging from 10 to 15 km, so if 
pure quark-stars exist, they should in principle be distinguishable from ordinary 
neutron stars by observation. Interestingly, several observations of neutron star 
candidates possessing properties that are in apparent contradiction with current 
neutron star models exist, most prominently the object RX J1856. 5-3754, whose 
inferred mass-radius relation (see Ref. [92, 93]) is also shown in Fig. 4.8. The 
exciting possibility that with the observation of RX J1856. 5-36754 one has a direct 
confirmation of the existence of quark-stars is still a matter of ongoing debate, 
as there are propositions claiming that the available data can be described by an 
ordinary neutron star; future observations are therefore needed to settle the issue 
of the exact nature of RX J1856. 5-36754. 

Finally, it should be kept in mind that the outermost layers of such a possible 
quark star consist of hadronic matter, giving rise to what is usually named a 
"hybrid star". The details of the star structure will depend sensitively on the 
hadronic equation of state [45], which unfortunately does not match on naturally 
to the EOS of the quark-gluon plasma phase calculated above [46]. Therefore, to 
make clear predictions for the phenomenology of compact stars one has to find 
a way to describe the intermediate range of densities in the equation of state, 
possibly by considering effective field theory models [97]. 

4.6 Application 2: Expansion of the quark-gluon 
plasma 

As a second application, I shall now describe the longitudinal expansion of the 
matter produced in a central collision of identical nuclei, using Bjorken's relativis- 
ts hydrodynamic model and the quasiparticle EOS for the quark-gluon plasma 
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calculated above. The starting point is the energy-momentum tensor of a rela- 
tivistic fluid with local energy density e(x) and pressure p(x), 

T^u = (e + p)U^U u - g^p, (4.13) 

where is the four- velocity of the fluid obeying U 2 — 1. Energy, momentum 
as well as the baryon number iV^ = nU^ 1 of the system should be conserved, 
therefore 

= 0, &>Np = 0. (4.14) 
Introducing proper time r and space-time rapidity r] defined as 

t = Vt 2 - z 2 , 7]=^-\n 1 —^- (4.15) 

and neglecting the effect of the transverse expansion as well as assuming that the 
properties of the fluid are invariant under Lorentz boosts in the z-direction this 
implies [98] 

e = e(r), p= P (r), T — T(t), = - (t, 0, 0, z) . (4.16) 

r 

With the use of the expressions 

d li T = U ll , d v U li = -(g IM ,-U IM U v ) (4.17) 
r 

and g^ u = (1, 0, 0, 1), the conservation equations can be brought into the form 

(£) = 0, d»N» = 0, (4.18) 

implying that the entropy density s per number density n is conserved by the 
longitudinal expansion. 

Using the quasiparticle model EOS one can calculate the curves of constant 
s/n in the T, /z-plane, shown in Fig. 4.9. Surprisingly, the curves are very well 
approximated by straight lines having slopes consistent with the tree-level ap- 
proximation 

rp3 ( 327T 2 , 2$TT 2 \ 

s_ 1 {— + _T 74tt 2 

n ~ ^T 2 ~ ^"9T' { ' 

as can be seen in Fig. 4.9; only at very large \i the tree level approximation 
seemingly starts to break down. 

Provided one has knowledge about the value of T, \i that the system reaches at 
freeze-out (e.g. by analyzing experimental data using thermal model descriptions 
of particle production [99]) as well as an estimate of the initial energy density 
created after heavy-ion collision, the above results allow an estimate of the initial 
temperature and chemical potential. For example, Braun-Munzinger et al. 
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Figure 4.9: Curves of constant s/n = 100, 50, 30, 10, 5 for the HTL model (dots). 
Also shown are the tree-level results (dashed straight lines) for the corresponding 
values of s/n. 



obtained (T, /x) ~ (174MeV, 15MeV) as freeze-out values for the central Au-Au 
collisions at RHIC with y/s = 130GeV, which corresponds to a value of so/n ~ 
100. Assuming furthermore an initial energy density of e = 20GeV/fm 3 ~ 175T C 4 
[99] and again T c = 172MeV one obtains (T ,/x ) = (344MeV, 27MeV) as initial 
temperature and chemical potential. Consequently, one can calculate the time 
behavior of the system by e.g. integrating Eq.(4.18), 

n = ™, (4.20) 
r 

where n = n(T(/i ), /i ) and r is expected to be on the order of r ~ lfm/c. In 
this model, one finds that the system freeze-out time is given by 17 ~ 12fm/c. It 
is clear however, that the above result is only a rough estimate of the time the 
system remains in a quark-gluon plasma phase since one has to take into account 
viscous effects of the fluid as well as the transversal expansion of the system that 
cannot be neglected at late times. 



4.7 Summary 

In this chapter, I have set up a possible extension of the HTL quasiparticle model 
that includes the full plasmon effect, however at the price of having to introduce 
another parameter. The values of the parameters used to fit the quasiparticle 
models to the Nf = 2 lattice data were given and it was shown that the resulting 
running coupling turns out to be comparable to the 2-loop perturbative coupling. 
I then investigated the quark-number susceptibilities at vanishing chemical po- 
tential and showed that they agree very well with independent lattice data; for 
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the HTL model, it was also shown that the higher-order derivatives represent 
the trend of another independent lattice study. Furthermore, it turned out that 
the quasiparticle pressure can be very well described by these [i = susceptibili- 
ties, even up to rather large values of chemical potential. Moreover, I calculated 
the lines of constant pressure starting from the critical temperature at \x = 0, 
which were shown to be compatible with lattice studies in the validity region 
of the latter. Extrapolating these results to large chemical potentials and small 
temperatures I obtained an equation of state for cold dense matter which should 
represent an improvement over existing perturbative results. As a first applica- 
tion, I then derived the mass-radius relationship of so-called quark-stars, finding 
results that fit nicely with the values inferred for RX J1856. 5-3754, a possible 
candidate for a quark star. As a second application, I calculated the expansion 
of a quark-gluon plasma created after a heavy-ion collision in the Bjorken model, 
finding that the system cools along nearly straight lines in the T, /i-diagram. Fi- 
nally, I estimated the time that the system remains in the quark-gluon plasma 
phase using freeze-out values of the temperature and chemical potential inferred 
for Au-Au collisions at RHIC. 



Chapter 5 

Results independent of the lattice 



In the last chapters I used an effective coupling to describe lattice data at van- 
ishing chemical potential which eventually gave a plausible equation of state for 
arbitrary chemical potentials. Although this approach seems to work very well 
it is somewhat unsatisfactory to be limited to available lattice data and a phe- 
nomenological ansatz for an effective coupling. Therefore, I want to investigate 
in this chapter how the results of the previous chapters are modified if one re- 
nounces all lattice input safe one number, namely the ratio T c /A^g, which for 
Nf = 2 I take to be given by 0.49, as was used before. Instead of the effective 
coupling I will use the standard two-loop running coupling of Eq.(1.6) with jl 
ranging between nT and 4nT at /x = 0, therefore testing the theoretical error 
of the model predictions. Furthermore, since the quasiparticle models represent 
an implementation of the thermodynamic quantities which is correct at least to 
leading perturbative order, one is then able to make predictions of these quan- 
tities "uncontaminated" by lattice artifacts using the machinery of the previous 
sections. I will then compare these predictions to both perturbational approaches 
and lattice results (whenever available) as well as those obtained in the previous 
chapters. 



5.1 Setup 

Using the HTL and NLA quasiparticle models for the pressure and the two-loop 
perturbative coupling, the expression for the entropy becomes that of Blaizot 
et.al [29]. Evaluating the entropy as a function of T by using T c /Aj^g = 0.49 
for Nf = 2 (which differs from the value used in the original reference) one can 
compare the result to lattice data for Nf = 2 from Ref. [18] (used for the fits in 
the previous sections) and Ref. [19] in Fig. 5.1a. As can be seen, the NLA model 
results are in general agreement with the lattice data (especially since a rigorous 
extrapolation of the lattice data to the continuum limit is still missing) while the 
HTL model seems to predict slightly too small values for the entropy. 
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Figure 5.1: Two-flavor entropy (a) and pressure (b) for HTL (long-dashed lines) 
and NLA quasiparticle models (full lines for c\ — 1, dotted lines for ca = 2, |), 
where the renormalization scale has been taken to be p = 4nT and p = nT, 
respectively. Also shown are Nf = 2 lattice results from Refs. [18] (triangles) and 
[19] (boxes), respectively. 

The pressure of the respective models is then given by Eq.(3.9), where B is 
fixed in such a way that the pressure at T = T c is equal to the pressure of a free 
gas of 3 massless flavors of pions, 

p(T c ,0)/T c 4 = ^ 0.329. (5.1) 

This completes the setup of the model. The resulting pressure at \i = as a 
function of the temperature is shown in Fig. 5.1b; although the model pressure is 
not as close to the lattice data as in the previous chapters where a phenomeno- 
logical fit of the coupling was used, it is still remarkable that the resulting band 
turns out to cover the lattice results and not much more while the lattice input 
has been minimized to one number. 



5.2 Results 

The characteristic curves again take shapes resembling ellipses as was the case 
in the previous chapters. Calculating the quark-number susceptibilities one finds 
the results shown in Fig. 5.2. As can be seen, the agreement with lattice data 
is still fairly good at high temperatures, whereas at temperatures comparable to 
T c the variations from varying p by a factor of 2 around 2nT get rather large. 
Note however that the results from the NLA models for the susceptibilities clearly 
differ from those obtained originally by BLAlzOTei al. [100]. The discrepancy is 
due to the fact that here the fermionic contributions were divided into soft and 
hard part similar to the gluonic contributions while in the original reference the 
fermions have been treated differently. As a consequence of treating the fermionic 
sector as in Ref. [79], the NLA model susceptibilities for Nf = 2 lie only slightly 



5.2. RESULTS 



47 



x/xo 




1.5 2 2.5 3 3.5 4 



T/T c 

Figure 5.2: Two- flavor susceptibilities for HTL (long-dashed lines) and NLA 
quasiparticle models (full lines for ca = 1, dotted lines for c\ = 4, |), where 
the renormalization scale has been taken to be p = 4nT and p = nT, respec- 
tively. Also shown are Nf = 2 lattice results from Ref. [85] (full triangles) and [40] 
(open boxes), respectively. The gray band corresponds to results from strict per- 
turbation theory to order a^lna [86] (see text for details). 



higher than the HTL results, as shown in Fig. 5.2. Therefore, the NLA results 
are in good agreement with existing lattice data [85, 40] while the original NLO 
results did not have any overlap with the data. I also show the results for the 
quark-number susceptibilities from strict perturbation theory to order In a s by 
Vuorinen [86, 87] in Fig. 5.2; the band shown corresponds to varying p = irT 
to p = 4nT while the unknown coefficient is taken close to zero [87]. 

When calculating the slope of the lines of constant pressure one obtains 
T c jgr ~ -0.055/ - 0.065 for the HTL quasiparticle model for p = txT (first 
value) and p = 4tiT (second value). The variations for the NLA model are 
slightly larger since one also has to consider the variation in the coefficient c\; 
one finds T c -^ ~ —0.052/ — 0.065 for the NLA model when varying p as above 
and ca by a factor of 2 around 1. Therefore, also the slopes of constant pressure 
turn out to be very similar to what was found in chapter 4, indicating that the 
agreement between these results and lattice calculations for Nf = 2 is a quite 
robust prediction of the quasiparticle models and in fact nearly independent of 
the fit data used in the previous chapters. 

Extrapolating the Nf = 2 quasiparticle results to the region of small temper- 
atures one can calculate an estimate for the critical density /x where p(0, /i ) = 0. 
Again varying p one finds for the HTL model /z — 3.49T C and /i = 2.93T C for 
p = iiT and p = 4tiT, respectively. For the NLA model the results turn out to 
be very similar, with the only difference being that the upper prediction of /Iq is 
slightly increased to /zq — 3.52T C . Comparing these to the results from the fitted 
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Figure 5.3: Two-loop (a) and one loop (b) fits (full line) to the quasiparticle 
model coupling at T = (points). 



quasiparticle models one finds that for ft = nT the estimates for /x are noticeably 
higher in the models not fitted to lattice data at \i = 0, while for larger values of 
ft the results become very similar (i.e. for ft = AnT the HTL model result is only 
about 2% lower than in chapter 4). 

Finally, it is also interesting to investigate closer what form the coupling 
constant takes in the region T -C £t- For example, one can try to describe 
a s (T = 0,/i) by the standard two-loop result with ft — yft, where y is to be 
determined by a least-square fit. This approach turns out to be not too successful, 
as can be seen in the exemplary plot of Fig. 5.3a, where such a fit has been tried 
on the HTL model result. Remarkably, the situation is very different when using 
the one-loop form for the running coupling Eq.(1.4) instead of the two-loop result. 
For the one-loop coupling, the fits with ft — yft typically look like the one shown 
in Fig. 5.3b, both for the HTL model as well as for the NLA models. Moreover, 
it turns out that the fit constant c is always proportional to the scale used at 
p — 0, i.e. when the renormalization scale was ft = AnT at ft = the scale that 
fits the coupling best at T = is given by ft ~ An\i 0.37 while for an original scale 
of ft = nT one finds ft ~ 7r/i0.33 at T = 0. 



5.2.1 Notes on changing T c /Ams 

As stressed before, the above results depend only on one external parameter, 
namely the value of T c /A^. It is therefore interesting to investigate how the 
results derived above will change when this parameter turns out to be smaller or 
bigger than 0.49. To this effect it is sufficient to see that the running coupling 
only depends on the scale through ft/Aj^ and accordingly 

Increasing the value of T c /Aj^g therefore amounts to increasing the ratio renor- 
malization scale over critical temperature. As an effect, the entropy and the sus- 
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ceptibilities will rise, whereas the slopes of constant pressure will become more 
negative. Finally, the estimate for the critical density will also take smaller values. 

5.3 Summary 

In this chapter I investigated the results from various quasiparticle models when 
using only one number, namely the value of T c /Aj^g as an input. It turns out that 
the pressure and the quark number susceptibilities at \i = are still described 
fairly well. The results for the slopes of constant pressure are close to what has 
been found in the previous chapters as long as p, is not much smaller than 2ttT, 
as is the case for the estimates of the critical density at T = 0. Furthermore, it 
turns out that the resulting coupling at T = can be very well described by the 
standard one-loop running coupling. I conclude by pointing out how the results 
derived in this chapter have to be modified when T c /A^ changes from the value 
adopted here. 
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Chapter 6 

The anisotropic quark-gluon plasma 



In the previous chapters I have limited myself to the consideration of isotropic 
systems. However, as pointed out in the introduction, the only system above 
the deconfinement transition one can (probably) study on earth, namely that of 
a fireball following an ultrarelativistic heavy-ion collision, is rather anisotropic 
in nature. Therefore, it is interesting to investigate which effects the presence 
of an anisotropy has on the dynamics of the system and which differences one 
can expect when comparing to the usually studied isotropic case. Note that I 
will consider only systems which are anisotropic in momentum space, therefore 
taking a snapshot of the system at some short time after the collision and before 
the system had time to convert the anisotropy from momentum space to config- 
uration space. Although eventually one would require a full treatment including 
anisotropies in configuration space, I will show that one can learn a lot about 
anisotropic systems by considering momentum anisotropies alone; moreover, a 
treatment of configuration space anisotropies would involve a full non-equilibrium 
quantum field theoretical description of the system, which despite recent progress 
is not available yet. 

In this chapter, I will therefore analyze the collective modes of high-temperature 
QCD in the case when there is an anisotropy in the momentum-space distri- 
bution function. I will present results (published in Ref. [101]) for a class of 
anisotropic distribution functions which can be obtained by stretching or squeez- 
ing an isotropic distribution function along one direction, thereby preserving a 
cylindrical symmetry in momentum space. 

6.1 Gluon self-energy in an anisotropic system 

Generalizing the classification of quasiparticles in an isotropic quark-gluon plasma 
in chapter 2 one first derives the HTL resummed gluon self-energy in an anisotropic 
system within semi-classical transport theory [58, 59, 60], which has been shown 
to be equivalent to the HTL diagrammatic approach [102]. Within semi-classical 
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transport theory partons are described by their phase-space densities and their 
time evolution is given by Vlasov-type transport equations [103, 104]. Concen- 

1/2 

trating on the physics at the soft scale, k ~ a s T <T, the magnitude of the field 
fluctuations at this scale is A ~ a l J 4 T and derivatives are of the scale d x ~ aJ 2 T . 
With this power-counting a systematic truncation of the terms contributing to 
the transport equations for soft momenta can be realized. 

At leading order in the coupling constant the color current J M induced by a soft 
gauge field A^ with four-momentum K = (u, k) can be obtained by performing a 
covariant gradient expansion of the quark and gluon Wigner functions in mean- 
field approximation. The result is 



where V 1 = (l,k/u;) is the velocity of the hard plasma constituents, Sn a (p, X) 
is the fluctuating part of the gluon density, and 8f+(p,X) and 5f"(p,X) are the 
fluctuating parts of the quark and anti-quark densities, respectively. Note that 
5n a transforms as a vector in the adjoint representation (5n = Sn a T a ) and 6f± 
transforms as a vector in the fundamental representation (5f± = 5f±t a ). 

Neglecting the collision terms which only enter at subleading order, the quark 
and gluon density matrices above satisfy the following transport equations: 



where Dx = dx +i(4:na s ) 1 ^ 2 A(X) is the covariant derivative and n(p) and /±(p) 
are the initial anisotropic gluon and quark distribution functions (which in the 
isotropic limit would correspond to the Bose-Einstein and Fermi-Dirac distribu- 
tion functions, respectively). 

Solving the transport equations (6.2) and (6.3) for the fluctuations Sn and 
Sf± gives the induced current via Eq.(6.1), 



where U(X, Y) is a gauge parallel transporter defined by the path-ordered integral 




[V-D x ,6f ± (p,X)\ 
[V-D x ,8n(p,X)} 



T{^a s f/ 2 V,F^{p,X)d v f ± { V ) , 
-{Aita s ) l l 2 V,F^( Pl X)d u n(3>) , 



(6.2) 
(6.3) 




x / dr U(X, X - Vr)F af3 (X - Vr)U(X - Vr, X) , (6.4) 




Y 



(6.5) 



F a p = d a Ap — dpA a — i(47ro; s ) 1 / 2 [A At , A v \ is the gluon field strength tensor, and 



h(p) = 2Nn(p) + N f (f + (p) + f.(p)) 



(6.6) 



6.1. GL UON SELF-ENERG Y IN AN ANISOTROPIC SYSTEM 



53 



is a particular combination of the gluon and quark distribution functions. Ne- 
glecting terms of subleading order in a s (implying U — > 1 and F a p — + d a Ap — 
dpA a ) and performing a Fourier transform of the induced current to momentum 
space one obtains 



J UK) = 4™ s / ^V% )h ( P ) (g, fi - ^^-) , (6.7) 



where e is a small parameter that has to be sent to zero in the end. 

From this expression of the induced current the self-energy is obtained through 
the relation 

w ^ - §M • < 6 - 8 > 

which gives the gluon self-energy of an anisotropic system, 

iV(K) = A,a s J |^ V&Mp) - K^vh) ■ ^ 

Note that this result can also be obtained using diagrammatic methods if one 
assumes that the distribution function is symmetric under p — > — p [102]. Fur- 
thermore, in the isotropic limit h(p) — > hi SO (p), one recovers the well-known HTL 
gluon self-energy from Eq.(6.9). 

After a little bit of algebra one finds that this tensor is symmetric, IP" (If) = 
n" M (if), and transverse, K^H^^K) = 0, if the distribution function vanishes on 
a two-sphere at infinity, linip^oo h(p) = 0, since e.g. 

K^ = -Ava s k 3 J -0^&h(p). (6.10) 

Therefore, only the spatial components of the self-energy are needed for the dis- 
persion relations of the quasi-gluons, which can be shown as follows: in the linear 
approximation the equations of motion for the gauge fields can be obtained by 
expressing the induced current in terms of the self-energy 

JUK) = W"{K)A"{K) , (6.11) 

and plugging this into Maxwell's equations 

-iK,F^{K) = JUK) + JZJK) , (6.12) 



to obtain 



[K 2 g^ - K»K V + W V {K)\A V {K) = -J^JK) , (6.13) 



where J" xt is an external current. Using the gauge covariance of the self-energy 
in the HTL-approximation one can write this in terms of a physical electric field 
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by specifying a particular gauge. Using the temporal axial gauge (A = 0) which 
has already been used in the isotropic case in chapter 2 one obtains 

[{k 2 - u 2 )S ij - k l W + H ij (K)]E j (K) = {&-\K)) ij E j {K) = ito J l eKt {K) . (6.14) 

Inverting the propagator allows one to determine the response of the system to 
the external source 

E\K)=iuW{K)jUK). (6.15) 

The dispersion relations for the collective modes can then be obtained as in 
chapter 2, by finding the poles of the propagator A^(K). 

6.1.1 Tensor decomposition 

Since an anisotropic system possesses a preferred spatial direction 1 , the gluon 
propagator cannot simply be decomposed into a transversal and longitudinal 
part as was the case for the isotropic system treated in chapter 2; therefore one 
needs to construct a new tensor basis. 

As mentioned above, the gluon self-energy is symmetric and transverse; as 
a result not all components of LP" are independent and one can restrict the 
considerations to the spatial part of IP", denoted IP 7 . One therefore needs to 
construct a basis for a symmetric 3-tensor that - apart from the momentum k % 
- also depends on a fixed anisotropy three- vector h\ with n 2 = 1. Following 
Ref. [66] one first defines the projection operator 

A ij = gij _ tftf/k 2 , (6.16) 

and uses it to construct n l = A^rV which obeys h ■ k = 0. With this one can 
construct the remaining three symmetric 3-tensors 

B" = k l k 3 jk 2 (6.17) 

C ij = nW/n 2 (6.18) 

D ij = k i h j + k j fi\ (6.19) 

The basis spanned by the four tensors A, B, C, and D therefore allows one to 
decompose any symmetric 3-tensor T into 

T = aA + bB + cC + dB; (6.20) 

furthermore, the inverse of any such tensor is then given as 

(a + c)B-a-\bc-h 2 k 2 d 2 )C-dT> 
= a A + b{a + c)-nW ' (6 - 21) 



x If I were to consider a system where also the cylindrical symmetry is lost there would be 
two preferred spatial directions and the tensor basis for this system would accordingly be more 
complicated, but still can be constructed by the method proposed here. 
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as a short amount of algebra reveals. With the above tensor decomposition one 
is now able to decompose the propagator into modes similar to what is done in 
the isotropic case. For this, one first needs to decompose the self-energy into 
"self-energy structure functions", which are the equivalents of the transversal and 
longitudinal parts of the self-energy in the isotropic case. 



6.1.2 Self-energy structure functions 

From Eq.(6.9) the spatial part of the general gluon self-energy tensor is found to 
be 

U ij (K) = -4na s [ -^-v l d l h(p) (V + —^—] . (6.22) 



(2tt) 3 ^ j V K-V + ie 

Since the collision term was set to zero, the distribution function h(p) is com- 
pletely arbitrary at this point, so in order to proceed one needs to assume a 
specific form for the distribution function. In what follows I will require that 
h(p) can be obtained from an (arbitrary) isotropic distribution function by the 
rescaling of one direction in momentum space, 

Hp) = Hp) = #(0 h iso (Vp 2 + £(p • fi) 2 ) , (6.23) 

where n is the direction of the anisotropy, £ > — 1 is a parameter reflecting the 
strength of the anisotropy and iV(£) is a normalization constant 2 . Note that 
£ > corresponds to a contraction of the distribution in the n direction (shown 
in Fig. 6.1) whereas —1 < £ < corresponds to a stretching of the distribution 
in the n direction. 

The normalization constant iV(£) has been introduced to preserve the relation 

d 3 n , . . , 

jw = /I(P) (6 - 24) 

between number density n and the momentum-space distribution function also 
in the case where the latter is anisotropic. N(£) is then simply determined by 
requiring n to be the same both for isotropic and anisotropic systems, 

/ Wf kM = I (2^ (P) = m I Wf kis ° (vV + to-fi) 2 ) (6-25) 
and can be evaluated to be 

N(0 = y/T+Z (6.26) 



2 In the original reference [101] this normalization constant has been set to one since it does 
not affect the collective modes qualitatively. However, the correct value has to be used in 
order to make quantitative predictions of physical observables (e.g. applying the results of this 
chapter to the calculation of the heavy quark energy loss in chapter 7) . 
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Figure 6.1: Contour plot of an anisotropic version Eq.(6.23) of the Fermi-Dirac 
distribution function with positive anisotropy parameter. The anisotropy vector 
is taken to be along the p^-direction. 



by performing a change of variables to p 

p 2 = P 2 (1 + £(v ■ n) 2 ) . (6.27) 

The same change of variables then also allows one to simplify (6.22) since then 
it is possible to integrate out the |p|-dependence, giving 

where 

m D = / dpp , 6.29 

7T J dp 

which corresponds to Eq.(2.4). Using the tensor basis from above, one can then 
decompose the self-energy into four structure functions 

U ij = aA 11 + (3B ij + -fC 13 + 5D ij , (6.30) 

which are determined by taking the contractions 

km ij k j = k- i . 

n { U ij k j = n 2 k 2 5 , 
fjll ij n j = n 2 (a + 7) , 
TrU ij = 2a + /5 + 7 . (6.31) 

Since the form of the structure functions is somewhat unwieldy, the integral 
expressions following from Eq.(6.31) have been relegated to appendix C. 




Figure 6.2: Real and imaginary parts of a/m 2 D as a function of real ui/k (a) and 
(b) the real part of a/m 2 D for uj/k = iT/k with 6 n = n/4 and £ = {0, 1, 10} (full 
line, dotted line and dashed line, respectively) in both cases. 

6.1.3 Evaluation of the structure functions 

The four structure functions a,/3,7,<5 depend on the Debye mass mo, the fre- 
quency and spatial momentum ui and k, the strength of the anisotropy £ and the 
angle between the spatial momentum and the anisotropy direction k • n = cos^. 
For symmetry reasons one can restrict < 9 n < | in the following. In the 
isotropic limit (£ — > 0) the structure functions a and f3 reduce to the isotropic 
HTL self-energies and 7 and 5 vanish, 

a(K,0) = U T (K), 
7(^,0) = 0, 

5{K,0) = 0, (6.32) 

where Ut(K) and IIl(A;) are given by Eq.(2.3). Note that for finite £ the analytic 
structure of a, ft, 7, 5 is the same as for Ilx,, rhr in the isotropic case, namely there 
is a cut in the complex 00 plane which can be chosen to run along the real u axis 
from —k < to < k. For real-valued 00 the structure functions are complex for all 
to < k (corresponding to the Landau damping regime of chapter 2) and real for 
uj > k while for imaginary-valued uj all four structure functions are real-valued. 
As an example, a plot of the structure function a for real and imaginary values 
of uj, £ — {0, 1, 10}, and 9 n = tt/4 is shown in Fig. 6.2. 

Using these structure functions it is then possible to write the inverse propa- 
gator A _1 (fT) in terms of the above tensor basis 

A-\K) = (k 2 -u 2 + a)A + ((3- uj 2 )B + 7 C + 5B . (6.33) 



Applying the inversion formula (6.21) one obtains a decomposition for the gluon 
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propagator into modes 

A(K) = A A A + (k 2 -uj 2 + a + 7) A G B + [(/? - lu 2 ) A g - A A ]C - 5A G B , (6.34) 
where 

A A \K) = k 2 -oo 2 + a, (6.35) 
A G \K) = (k 2 -uj 2 + a + 1 )([3-uj 2 )-k 2 fi 2 5 2 . (6.36) 

In the following it is convenient to use the form 

A(K) = A A [A - C] + A G \{k 2 - u 2 + a + 7 )B + ((3 - u 2 )C - 5T>] , (6.37) 

so that the poles of the propagator correspond to the poles of and A<3, 
respectively. 

6.1.4 The quark propagator in an anisotropic system 

Using the equivalence of semi-classical kinetic theory and the HTL diagrammatic 
approach one can derive the quark propagator in an anisotropic system similar 
to the gluon propagator above. However, unlike the gluon propagator the quark 
propagator in an anisotropic system can be shown to behave very similar to its 
isotropic counterpart [105]. Consequently, the quark collective modes will look 
very similar to those derived in chapter 2 and therefore the quark propagator in 
an anisotropic system will not be treated here explicitly. 

6.2 Collective modes of an anisotropic quark- gluon 
plasma 

The dispersion relations for the gluonic modes in an anisotropic quark-gluon 
plasma are determined by setting A^ 1 = and A^ 1 = and then solving these 
equations using Eqs. (6.35,6.36). A subsequent comparison to the results of chap- 
ter 2 will show the qualitative differences with respect to the isotropic case. 

6.2.1 Static limit 

As a first test of how the momentum-space anisotropy in the distribution functions 
affects the properties of the system one can consider the response of the system 
to static electric and magnetic fluctuations. For this test, one needs to examine 
the limit uj — > of the propagators (6.35) and (6.36): approaching along the real 
us axis one finds that to leading order a ~ 7 ~ 0(u°), (3 ~ 0(u 2 ), and 5 ~ 0(iu); 
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therefore, one can define the four mass scales 

m 2 a = lim a , 



w^O {J 1 

2 



ml = lim -f3 , 



m„ = lim 7 , 

m\ = lim lm5 , (6.38) 

which will be used in the following. The static limit of the propagators A^ 
Eq.(6.35) and Ac Eq.(6.36) in terms of these masses is then given by the expres- 
sions 



A A L = k 2 + m 2 a (6.39) 
A G X = -^[(e + ml + m^ + m^-mj] . (6.40) 



Furthermore, Aq 1 can be factorized into 



A G X = -^(k 2 + ml)(k 2 + m 2 _) , (6.41) 



where 



and 



2ml = M 2 ± JM 4 - 4{m 2 Jm 2 a + m*) - mj) , (6.42) 



M 2 = m 2 a + ml + m 2 . (6.43) 

In the isotropic limit £ — > 0, most of the mass parameters vanish, m 2 a = m 2 = 
m 2 = m 2 _ — > 0, while m\ becomes equal to the Debye mass squared, m\ — > m 2 D , 
so that one recovers the results obtained in chapter 2. For finite £ it is possible 
to evaluate all four masses defined above analytically (the results for m a and mp 
are listed in Appendix C), but since the resulting expressions are quite unwieldy 
it is more straightforward to use a numerical evaluation. As an example, the 
angular dependence of m 2 a , m 2 + , and m 2 _ at fixed £ = 10 and £ = —0.9 is plotted 
in Fig. 6.3. In the case £ > (Fig. 6.3a), one can see that the scale m 2 + is bigger 
than m 2 D for small 6 n while being smaller than m 2 D for 6 n near n/2; indeed one 
finds that this corresponds to the screening of electrostatic modes in the isotropic 
limit, the only difference being an angular dependence of the screening length for 
nonzero £. 

Non- vanishing scales m 2 a and m 2 _, however, should naively correspond to a 
screening of magnetostatic modes, which is absent in isotropic systems. Indeed, 
the fact that anisotropic systems allow for non- vanishing magnetostatic screening 
masses was already found by Cooper et al. [106]; however, Cooper et al. 
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(a) (b) 




Figure 6.3: Angular dependence of m 2 a , m 2 + , and m 2 _ (dotted line, full line and 
dashed line, respectively) at fixed (a) £ = 10 and (b) £ = —0.9. 

neglected to consider the fact that the resulting mass squares might be negative 
(as is the case for small 6 n in Fig. 6.3a) and would therefore not correspond to a 
screening of the magnetic interaction. 

More precisely, the fact that these quantities are negative indicates that for 
£ > the system possesses an instability to transverse and "mixed" external 
perturbations associated with ra„ and m/L, respectively, as has now been veri- 
fied by BlRSE et al. [107] and discussed in a more general context by Arnold, 
Lenaghan and Moore [64]. As can be seen in Fig. 6.3a the transverse insta- 
bility is present for any 9 n ^ tt/2 while the mixed instability is only present for 
n < 0f xcd with 9f xcd depending on the value of f . 

For the case £ < 0, which is shown in Fig. 6.3b, the ^-dependence of the 
scales m 2 + and m 2 _ is exactly reversed with respect to the case £ > 0, as one 
would expect naively from Eq.(6.23). The scale m 2 a , however, which was strictly 
negative for £ > 0, now turns out to be strictly positive in the case of £ < 0, so 
the transverse instability vanishes in the case of £ < 0. For 9 n > tt/4 the mixed 
instability is still present, since the scale m_ is again negative. 

The role of these unstable modes will be treated later on in more detail. 

6.2.2 Stable modes 

In the non-static case, a factorization of A^ 1 similar to Eq.(6.41) can be achieved, 
which allows the determination of the dispersion relations for all of the collective 
modes in the system. 

Considering first the stable collective modes which correspond to poles of the 
propagator at real-valued u > k one factorizes A^ 1 as 

A G 1 = (uj 2 -n 2 + )(Lo 2 -nl), (6.44) 



where 

2Q 2 ± = Q 2 ± sJ& A - 4((a + 7 + k 2 )/3 - k 2 fi 2 5 2 ) , (6.45) 
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Figure 6.4: Angular dependence of E a (dotted line), E + (full line), and E- 
(dashed line) for £ = 10, and 9 n = {0, 7r/4, 7t/2} (plots (a), (b) and (c), respec- 
tively). 

and 

Ci 2 = a + (3 + -f + k 2 . (6.46) 

Since the quantity under the square root in (6.45) can be written as (a — (3 + 7 + 
k 2 ) 2 + 4k 2 n 2 5 2 (which is always positive for real uj > k), there are at most two 
stable modes coming from A^, while the remaining stable collective mode comes 
from the zero of A^ 1 . 

The dispersion relations for all of the collective modes are then given by the 
solutions to 

El = n 2 ± (E ± ), (6.47) 
E 2 a = k 2 + a(E a ) . (6.48) 

In the isotropic limit (6.32) one finds the correspondence E a = E + = Et and 
E- = El. Accordingly, for finite £, there are three stable quasiparticle modes 
with dispersion relations that depend on the angle of propagation with respect to 
the anisotropy vector, 9 n . The resulting dispersion relations for all three modes 
for the case £ = 10, and 9 n = {0, 7r/4, 7r/2} are plotted in Fig. 6.4. 

6.2.3 Unstable modes 

For non-zero £ the propagator also has poles along the imaginary uj axis 3 , that 
correspond to exponentially growing or decaying field amplitudes, depending on 
the sign of the imaginary part. The exponentially growing solutions correspond 
to the unstable modes of the system with u — > iT and T the real- valued solution 

3 Checking for poles at complex uj can be done numerically but no poles on the physical sheet 
have been found; indeed, for certain special cases one can prove that there are no other modes 
than the ones listed above [108]. However, there exist solutions on the unphysical sheet which 
come very close to the border of the physical sheet for sufficiently large values of the anisotropy 
parameter £ (see Ref.[108]). 
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Figure 6.5: T a (k) (dotted line) and r_(fc) as a function of k with (a) £ = 10 and 
6 n = vr/8 and (b) £ = -0.9 and 6 n = ir/2. 

to the equations 

a g x = (r 2 + ^)(r 2 + fi 2 ) = o, 

A^ 1 = (r 2 + k 2 + a) = 0. (6.49) 

It turns out that in contrast to the stable modes there is at most one solution 
for Aq^w = iT) = since numerically one finds that Q 2 > for all T > 0, while 
A^u; = iT) = has only solutions for £ > 0. Therefore, the system possesses 
one or two unstable modes depending on the sign of the anisotropy parameter. 

In Fig. 6.5a a plot of the growth rates T a (k) and r_(fc) with £ = 10 and 
Q n = is shown, reflecting the presence of two unstable modes; in Fig. 6.5b 
r_(fc) is shown for £ = —0.9 and 6 n = n/2. 



6.3 Small £ limit 

In the small-£ limit it is possible to obtain analytic expressions for all of the 
structure functions order-by-order in £, which provide a well-defined test case to 
the otherwise numerically conducted analysis. To linear order in £ one finds 



7 
6 



z 2 1 

— (3 + 5 cos 29 n )m 2 D - -(1 + cos 29 n )m 2 D 



z~ 2 (3 = -n L ( z ) + z 



+ ^n T (z) ((3 + 3cos2# n ) -2 2 (3 + 5cos2#„)) 
-(1 + 3 cos 26 n )m 2 D - (cos 26 n (2 - 3^ 2 ) + 1 - z 2 ) 



^(m T (z) -m 2 D )(z 2 - 1) sin 2 9 n: 



3A; 



(Az 2 m 2 D + 3U T (z)(l - Az 2 )) cos6 r , 



(6.50) 
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where z = ui/k and again k • n = cos# n . 



6.3.1 Static Limit 

Using the linear expansions from above and the well-known static limits TIl — ■> 
m 2 D and 11^ — > —muj/^Ak) one obtains analytic expressions for the mass scales 
(6.38), 

K = _|(l + cos20 n ) , 



ml = l + |(3cos2fl n + 2) , 



6 



~ 2 S • 2 n 

m 7 = - sm n , 

= —£—sm9 n cos9 n , (6.51) 

where m 2 = m 2 /m 2 D . Using these one obtains expansions to linear order in £ for 
m\ defined in (6.42) 

m\ = l + ^(3cos2# n + 2) , 

m 2 . = — ^-cos26 n . (6.52) 
o 

Clearly, the mass scales in the small-£ limit have the same qualitative behavior 
as their exact counterparts, as one can see by comparing the above expressions 
to Fig. 6.3. Therefore, one gains confidence that the small-£ limit conserves the 
basic features of the full anisotropic problem while having the advantage of being 
analytically tractable. 



6.3.2 Collective modes 

The collective modes in the small-£ limit are most easily found by an expansion 
of Eq.(6.34), 

A^ 1 = k 2 - u 2 + a = 

Aq 1 = (k 2 - uo 2 + a + -i){(3 - uo 2 ) = , (6.53) 

where a, f3, and 7 are given by (6.50) and S 2 can be ignored in the expansion 
since it is of order 0(£ 2 ). Solving the above equations one finds that again both 
the stable and unstable modes in the small-£ limit correspond qualitatively to 
the full problem. Note that there is again only one unstable mode coming from 
Aq 1 since (3{iV) > for all T > 0. 
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6.4 Large £ limit 

Another case where one is able to calculate analytic results for the structure 
functions is the limit £ — > oo. By considering the action of the distribution 
Eq.(6.23) on some smooth test function 3>(p), 



r 



lim 



(27T) 



r^(p)$(p) 



and choosing p • n to be parallel to the p 2 -axis one can rescale cos 
that the limit £ — > oo can be performed, 



r 



p 2 dpd<p 
(2tt) 3 



(pco 
psii 




(6.54) 
(6.55) 



(6.56) 



This means that in the limit £ — > oo Eq.(6.23) takes the form 

lim ft^(p) — > 5(p ■ n) y /i iso 

which corresponds to the extreme anisotropic case considered by Arnold, Lenaghan 
and Moore [64]. As a consequence, one can make use of this form by partially 
integrating Eq.(6.22) 



U ij {K) = Aiia s 



d 3 p h^(p) 

(27T) 



P 



5 lJ 



k l v j + k j v l (—uj 2 + k 2 )v l v j 



(6.57) 



-K-V-ie (-K-V-ie) 2 

and applying the techniques from Ref. [64] to obtain analytic expressions for the 
structure functions in the large £ limit. Using 

d 3 p hz(p) 



lim Ana s 



7T 



(2tt) 3 p 



m 



D~ 



one immediately finds 



U ij (K) 



m D 7i 



2n 



2ir 



k l v j + k j v l (—uj 2 + k 2 )v l v j 



+ 



-K-V-ie (-K -V -it) 



(6.58) 



(6.59) 



The remaining average over angles can easily be done, obtaining the structure 
functions through the contractions Eq.(6.31), giving 



a 



m D ix 



2/i 00 
-cot V n + 



sin 2 9 r , 



00 + 



l-LU 2 



yjoj + sm9 n ^u — sin 9 n 



(3 = 



-1 +oj— 



Co 2 + cos 29 n 



7 



mf,7r 1 - uj 2 
4 4sin 2 n 



(^ + sin0 n ) 3 / 2 (^-sin0 n ) 3 / 2 _ 

3 - 6£ 2 - 2(1 + uj 2 ) cos20 n - cos40 r , 



mjjir cos9 n 

o = 



4k sin 2 61 



UJ 



6 + 2cos2fl 
1 



a; 



(£ + sin0 n ) 3 / 2 (£-sin0 n ) 3 / 2 
2 ^ 2 + (l-2^ 2 )cos 2 n " 



(^ + sinfl n ) 3 / 2 (a>-sin0 n ) 3 / 2 



(6.60) 
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Figure 6.6: Angular dependence of the mass-scales m 2 a1 m 2 + ,m 2 _ (dotted line, full 
line and dashed line, respectively) in £ — > oo limit. 

It is then straightforward to extract the static limit of these structure functions, 
giving the mass scales m 2 a) m 2 + ,rn?i in the large-£ limit. The result is plotted in 
Fig. 6.6 as a function of the angle 6 n . As can be seen, the mass scales diverge 
at 6 n = 0, whereas m 2 _ and m 2 + coincidence at 6 n = tt/2 and m 2 a vanishes. Also, 
one recovers the results from Arnold, Lenaghan and Moore [64] that the 
scale m 2 _ is only negative for sin 9 n < 1/ V2 (corresponding to the presence of the 
electric instability in this reference) and the various growth rates of the unstable 
modes. Finally, one can show that the expressions of the structure functions for 
general £ converge towards the above analytic results when £ becomes large. 

6.5 £ — > — 1 limit 

To get some insight in the case of £ — > —1 where the distribution function 
Eq.(6.23) becomes extremely elongated one can again consider its action on a 
test function $(p), 

F = ,l m J (2^ (P)$(P) - (6 - 61) 
Performing a simple scaling p z — > -^== one obtains the interesting result 

F = Bm J dp, h iso ^Pl+P 2 )) $ ( ) (6.62) 

indicating that the argument of the test function <3> is pushed to very large values, 
where $ is usually taken to be vanishing. However, one can get a bit more insight 
by compactifying momentum space to a sphere and considering the coordinate 
map 

P - ^, (6.63) 
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with the volume element transforming to 

£p - ^. (6.64) 

One then considers the action of the distribution function on a test function $' 
in the primed coordinate system, 

F - J (^!) *'(P'). (6.65) 

Taking again p' ■ n = cos 9' one scales 

p' — > p y/l + £cos 2 9' (6.66) 

to obtain 

Since the isotropic distribution function consisting of Bose-Einstein and Fermi- 
Dirac statistical factors lifts the possible singularity at p' = 0, the only "danger- 
ous" term in the integrand is 1 + £ cos 2 9 which becomes zero for £ — > —1 and 
small angles 9'. Substituting therefore cos 2 9 — 1 — (l+£)x one can perform the 
limit £ — > — 1, obtaining 

(6.68) 

Therefore, in the limit £ — > — 1, Eq.(6.23) simply is given by 

/•oo 

lim /i ? (p)rf 3 p-*5 3 (p')rfV47r / dpp 2 h iso (p), (6.69) 

which upon insertion into Eq.(6.57) shows that in this limit the self-energy and 
consequently all structure functions vanish. Therefore, the limit £ — > — 1 un- 
fortunately does not correspond to the interesting physical situation of a line 
momentum distribution but rather to the free limit. 

6.6 Discussion of instabilities 

The existence of unstable modes in an anisotropic quark-gluon plasma (first stud- 
ied by Mrowczynski in a series of papers [58, 59, 60]) may play an important 
role in the dynamical evolution of the quark-gluon plasma, since their exponential 
growth can lead to a more rapid thermalization and isotropization of the plasma. 



6.7. SUMMARY 
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Recently, Mrowczynski and Randrup [63] have performed a phenomenolog- 
ical estimate of the growth rate of the instabilities for scenarios relevant to the 
quark- gluon plasma produced at RHIC or LHC, using however a different im- 
plementation of the anisotropy than Eq.(6.23). They found that the degree of 
amplification of the instabilities is not expected to dominate the dynamics of the 
quark-gluon plasma, but instead their effect would be comparable to the contri- 
bution from elastic Boltzmann collisions. However, they also pointed out that 
if a large number of unstable modes would be excited then their combined ef- 
fect on the overall dynamics could well be significant. More recently, Arnold, 
Lenaghan and Moore [64] investigated the case corresponding to the £ — > oo 
limit from above, arguing that it drastically modifies the scenario of "bottom-up 
thermalization" advocated by Baier, Mueller, Schiff and Son [53], which 
would then have to be replaced by a different scheme. 

Furthermore, the presence of instabilities will in general prohibit the calcu- 
lation of physical quantities in a perturbative framework, since they correspond 
to unregularized singularities of the propagator. This problem as well as special 
exception to it will be discussed in more detail in chapter 7. 

Because of these findings, it is probably fair to say that a more intensive 
study of the instabilities of anisotropic systems and especially their saturation 
and effect on the thermalization will be of great interest in the near future. 

6.7 Summary 

In this chapter I calculated the gluon self-energy in a system that has an anisotropy 
in momentum space. I presented a tensor basis that can be used to write the 
gluon self-energy as the sum of four "structure functions" and showed that in the 
isotropic limit two of these structure functions vanish while the others correspond 
to the longitudinal and transversal part of the self-energy of chapter 2. Consid- 
ering the static limit of the self-energy it was found that for anisotropic systems 
there is in general a non-vanishing magnetostatic mass-scale, but since the as- 
sociated mass squared is negative, this corresponds to the existence of unstable 
modes. 

Furthermore, solving the dispersion relations three stable and either one or 
two unstable modes were found, depending on whether the distribution function 
is stretched or squeezed. I also considered the limit of weak anisotropies and two 
cases where the anisotropy gets extremely strong, deriving analytic expressions 
for the structure functions which serve as verification of the general result and 
provide the basis for further analytic studies. I conclude by reviewing the current 
knowledge about the effects of the unstable modes on the thermalization of the 
quark-gluon plasma produced in heavy-ion collisions. 
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Chapter 7 

Energy loss of a heavy parton 



In this chapter I will make use of the resummed self-energies derived for isotropic 
and anisotropic systems obtained in chapters 2 and 6 to calculate the collisional 
energy loss of a heavy parton in a deconfined QCD plasma. In the context of 
heavy-ion collisions, a theoretical understanding of the heavy parton energy loss 
in a quark-gluon plasma is important for a proper comparison with experimental 
results, e.g. for the inclusive electron spectrum (measured recently at RHIC 
[109]) and for jet suppression rates. In the following I will focus on the collisional 
energy loss (for a review on the radiative energy loss see e.g. [110, 111]) using the 
technique of Braaten and Thoma [112, 113] which gives the complete leading- 
order result for this quantity. 

The main idea behind this technique is to consider independently the contri- 
butions from soft (involving momenta q ~ m D ) and hard (q ~ T) gluon exchange, 
which are separated by some arbitrary momentum scale q* that cuts off the UV 
and IR divergences of the soft/hard contributions, respectively. Moreover, it was 
found that in the weak-coupling limit a s C 1 the condition <C q* <T could 
be used to expand the resulting integral expressions for the soft and hard contri- 
butions further, giving an analytic result for the energy loss independent of the 
separation scale q*. However, in the case of QCD where the coupling becomes 
quite large, this approach may give unphysical results for energy loss, as will be 
shown later on. 

It turns out that when not expanding the integral expressions with respect 
to the condition m D <C q* C T, in the isotropic case both the hard and soft 
contributions to the energy loss always stay positive, even for very large coupling, 
as has been shown in QED [114]. However, this comes at the expense of giving up 
independence of the complete result on the scale q*, which then has to be fixed 
somehow. Fortunately, it turns out that in the weak coupling limit the scale 
dependence becomes very small, too, so that unless q* is taken to be very large 
(q* ^> T) or very small (q* <C mo) one recovers the result found by BRAATEN and 
THOMA. When the coupling is increased one can then still fix q* by the "principle 
of minimal sensitivity", which means that q* is chosen such that the energy loss 
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(and therefore also the variation with respect to q*) is minimized. As will be 
shown, this procedure eliminates the unphysical result of negative energy loss in 
the isotropic case and also provides an estimate on the theoretical uncertainty of 
the result by varying q* by a fixed amount around the minimizing value. 

Moreover, it will be shown that the technique by Braaten and Thoma 
can also be applied to anisotropic systems. This is in fact non-trivial since the 
presence of plasma instabilities in general could render the calculation of the 
soft part divergent; however, the collisional energy loss turns out to be protected 
by an interesting mechanism dubbed "dynamical shielding" [114], rendering the 
calculation safe at least to leading order, as will be shown in the following. 

7.1 Soft contribution and effect of instabilities 

The soft contribution to the collisional energy loss can be calculated using classical 
field theory methods. Starting from the classical expression for the parton energy 
loss per unit of time given by 

ft = Re / ^ X J extP0 • KAX) , (7.1) 

where a is a color index, X = (t, x), and J ex t is the current induced by a test 
parton propagating with velocity v: 

J a cxt (X) = C°v<5< 3 >(x-vt), 

Jcxt(Q) = (27r)Q<V%-q-v), (7.2) 
with Q = (cu, q) and Q a being the color charge. Using 

2&(Q) = *" " A?(Q)) W) , (7.3) 

where Aq denotes the free propagator, and Fourier transforming to coordinate 
space one obtains 

E^(X) = i Q a J (q . V )(A«(Q) - A«(Q)y e^-h-W , (7.4) 

where Q = (q- v, q). Using the above equation one finds for the soft contribution 
Eq.(7.1) to the heavy parton energy-loss 

- ^ = Q 2 Im J (q . v y [A«(Q) - A?-(Q)] ^ . (7.5) 
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7.1.1 Dynamical shielding of instabilities 

In the case of an anisotropic system unstable modes are present that manifest 
themselves by poles of the propagator in the static limit lim^o A" 1 = 0, as has 
been discussed in chapter 6. To simplify the argument, one can use the tensor 
decomposition of the propagator in chapter 6 and concentrate on the poles of A^ 
given in Eq.(6.34) only, while the poles of Ac can be treated in a similar way. 
Schematically the soft contribution to the energy loss Eq.(7.5) is then given as 

Since the structure function a for £ > is negative-valued, 

lim a(u, q) = M 2 (-l + iDCj) + 0(cu 2 ) , (7.7) 

in the static limit the integrand in Eq.(7.6) seems to contain an unregulated 
singularity at q = M, which therefore renders the result divergent. Thus the 
presence of instabilities in a system is believed to generically prohibit the cal- 
culation of quantities in a perturbative framework since those quantities will be 
plagued by unregulated divergencies [115]. 

The schematic integrand for the energy loss in Eq.(7.6), however, turns out 
to be "shielded" from this singularity, as can be seen by taking the limit 

J™ ^Sf _ q 2£2 + m 2 (-1 + iDu) ~* M 2 (u + iD) ~MW' ^'^ 

This limit exists as long as the coefficient D in Eq.(7.7) is non- vanishing, being 
somewhat similar to dynamical screening of the magnetic sector of finite temper- 
ature QCD. Therefore, the integral Eq.(7.6) may safely be performed, as long as 
d / 0. 

Proof of dynamical shielding I: small £ 

Since dynamical shielding relies crucially on the fact that the imaginary part of 
the structure function is of 0(u) when the real part is negative- valued, a simple 
proof of this assumption valid for the weak-anisotropy or small-£ regime will be 
given here. For A^ the proof is straightforward since Im A^ 1 = Im a and from 
Eqs.(6.50) one finds 



lim a = m 2 n 



-|(1 + cos20 n ) - l ^cj(l + M(i + C0S 2^ n ; 



(7.9) 



From this expression it is clear that the imaginary value cannot change sign for 
any positive value of £, whereas for negative value of £, the static limit of a is 
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positive- valued, so no instability exists in A^. To prove that Aq is also safe, one 
first uses its small-£ decomposition Eq.(6.53) to see that only the a + 7 term can 
become negative- valued. Using once again the expressions in Eq.(6.50) one finds 
after a little bit of algebra 



lim a + 7 = m D 

c2j — >0 



ITT 



-(-2 + 4 cos 2 n ) - — Cj 1 + ±(-4 + 10 cos 2 
6 4 



. (7.10) 



For £ > the instability is present for cos 2 # n > |, while the imaginary part 
may change its sign only for cos 2 ^ < ^; for £ < 0, the instability occurs for 
cos 2 n < \ while even for the "worst case" £ = — 1 the imaginary part can only 
change its sign for cos 2 ^ > ^5, so also the Ac contribution to the energy loss 
is protected by dynamical shielding. This completes the proof for the case of 
small £. 



Proof of dynamical shielding II: large £ 



Also for very large £ it can be shown analytically that the energy loss calculation 
is protected by dynamical shielding. Considering the A^ contribution first one 
can make use of the analytic expressions in Eq.(6.60) to find 



lim a 



m 2 D 7i 



-cot % 



sin J Or, 



(7.11) 



where again it can immediately be seen that the 0(u) contribution to the imag- 
inary part never vanishes. For Aq it turns out that a factorization similar to the 
small-£ case is possible, 



A G (oj,q)=u; 2 (q 2 + q 2 + )(q 2 + q 2 _), 



(7.12) 



where 



q 2 - = 



m D ir 



16 sin 2 

m 2 D 7r 
16 sin 2 

4 



2(1 + cos 2 n ) + 2a/ cos 4 #„ + 4cos 2 n ] + 0(u) 

4-2 cos 2 n - 2a/cos 4 n + 4 cos 2 n + 

, a — - (8 + 2 cos 2 n - 3A/2 a/(8 + 2 cos 2 n ) cos 2 n ) 
(4 + cos 2 n ) sm n V / 



100 



(7.13) 



As can be seen, g 2 > in the whole region < n < f , so one can limit the con- 
siderations to q 2 _ . There a little bit of algebra shows that q 2 _ < for cos 2 6> Tt > |, 
while on the other hand the imaginary part changes sign at cos 2 n — §, so all an- 
gles n 7^ f are protected through dynamical shielding, while for the angle n — \ 
the singularity at q = of Aq is rendered harmless by the volume element of the 
integration. This completes the proof for the case of large £. 
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7.1.2 Application of dynamical shielding 

Unfortunately a proof that shows that dynamical shielding protects the soft con- 
tribution to the energy loss of naked singularities is rather difficult for general 
£. Numerical evidence however suggests that dynamical shielding also works 
for general £, which is not covered by the above proofs. Therefore, one can be 
reasonably confident that the singularities coming from the unstable modes are 
rendered safe everywhere and one may proceed to evaluate the soft contribution 
to the energy loss. 

Using the tensor basis derived in chapter 6 to invert the propagator in Eq.(7.5) 
and taking the contractions 



v l A l3 v 3 
v l B ij v j 
v l C ij v j 
v l D ij v j 



= ^"(q-v) 2 , 

= (q-) 2 , 

= (n-v) 2 /n 2 , 

= 2(q-v)(n-v) 



(7.14) 



one obtains 



dx J 



soft 



Q 2 



Im 



d 3 q 

(2^)3 



u A A (Q) - 



or 



+ ujA g (Q) 
1 

+ 



CO 



UJ 



— (q 2 -uJ 2 + a + -f) + ({3- UJ 2 ) 



:a-y) 2 

n 2 



n ■ v 



h 2 



uj(q 2 



\21 



n 2 



25uo{h ■ v) 

(7.15) 



where dW/dx = v 1 dW/dt and uj = q • v. After some algebraic transformations 
involving scaling out the momentum one obtains 



dW \ 
~dx~) 



soft 



d 3 q 



Q 2 

T 

V m ./ (2tt) 3 q{l-UJ 2 ) 



-a 



(q 2 — q 2 uj 2 



uj 2 



a 



q 2 A + B 



q 4 C + q 2 V + S 



n • v 



h 2 



(7.16) 



where 



A 

B 

C 
V 
£ 



n • v 



{a + -f)-2uj(l-uj 2 )(h-v)5 , 

\2 



a + -f)(3-n 2 5 2 j (1 -to 2 - 

-uj 2 (1-uj 2 ), 

-c^ 2 (a + 7) + (l-w 2 )/5 

a + >y)[3 -h 2 5 2 , 



n • v 



-), 



(7.17) 
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Figure 7.1: Sketch of the directional dependence of the energy loss: n is the 
direction of the anisotropy, v the velocity of the parton and 9 the angle between 
these. 



with Cj = q • v and 5 = q5. Since all the momentum dependence has been made 
explicit and possible singularities coming from the instabilities are known to be 
shielded one can perform the q integration, obtaining 



dW \ Q 2 

dx ) soft v 



Im 



f dtt g 



Cj 2 ) 



—a 



2(1 -Cj 2 ) 
+F(q*)-F(0) 



a 



(7.18) 



where 



_/ x A, , Anln ± ^ 2 c \ \ AD-2BC . VV 2 -4CS + V + 2Cq 2 

F(q) = — In {-AC (Cq 4 + Vq 2 + £))+ In - — , 

V ; 4C v vy ;; 4CVV 2 - 4CS VV 2 -4CS-V- 2Cq 2 

(7.19) 

and a UV momentum cutoff q* has been introduced on the q integration. 



7.1.3 Behavior of the soft part 

To get some idea on how the soft part behaves as a function of the parameters 
it is useful to consider once again the limit of small £ in Eq.(7.18). It turns out 
that in this limit one obtains 



dW 



dx 



soft, small— £ 



dW 

dx 



soft,iso 



dW 

dx 



(vn) 



soft,£i 



dW 

dx 



soft, § 2 



(7.20) 



with 
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Figure 7.2: Soft contribution to the 
energy loss (full line) as a function of 
q*/m D for v = 0.5 compared to the 
Braaten-Thoma result (dashed line). 
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Figure 7.3: The contribution 
— {dW/dx) soft £ 2 as a function of 
q*/m D for v = 0.5. 



'dW\ Q\ f v duo T v 2 -cu 2 „ , (i-^ 2 y 2 + n, 

■uj 



Q T f 



da; /sofUso *> 2 V-„(27r) 



n r In 



2(1 -^ 2 ) 2 * n 

in 



2 -n, 



(7.21) 



being the isotropic result for general q* (calculated originally by Thoma and 
Gyulassy [116]). Note that Eq.(7.21) corresponds to the Braaten-Thoma result 
when expanding the logarithms under the assumption q* ^> m D . The respective 
results are compared in Fig. 7.2: while the results are identical for large q* /mn 
the Braaten-Thoma result becomes negative for small q*/m D whereas the unex- 
panded result obtained through numerical integration is positive for all values of 
q*/m D . 

The effects of the anisotropy are encoded in the functions (dW^/da;) soft ^ and 
(dW/dx) soit £ 2 which resemble that of (dW/dx) so{t iso , but since they consist of 
many more terms than the isotropic contribution they will not be listed here 
explicitly. The direction of the heavy parton enters the small-£ result only through 
the explicit term (v • h) 2 /v 2 . The function — {dW / dx) soft ^ 2 which multiplies this 
term and therefore controls the directional dependence of the soft part is plotted 
in Fig. 7.3. As can be seen in this figure, the function is negative for small q* /m D 
but becomes positive for large q*/m D . Since the function will eventually be 
evaluated at some q*/rriD that minimizes the overall energy loss depending on the 
velocity and the coupling constant, the directional trend of the soft part reverses 
once this value of q* /mn crosses the point where (dVF/dx) soft ^ 2 = 0. Therefore, 
one expects the energy loss to be peaked along n for some couplings and velocities, 
while being peaked transverse to n for others (see sketch in Fig. 7.1). 

For general £ it turns out that the results obtained in the small £ limit still 
hold qualitatively although quantitatively the predictions differ considerably once 
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£ becomes large. Fortunately, in the limit £ — > oo one can again use the analytic 
results for the structure functions of chapter 6 as a verification of quantitative 
results. 



7.2 Hard Contribution 



The hard contribution can be separated into two parts: one contribution coming 
from the scattering of the heavy parton on quarks in the plasma and another one 
that takes into account the scattering on gluons (corresponding to the tree-level 
diagrams shown in Fig. 7. 4). Assuming the velocity of the parton to be much 
higher then the ratio of the plasma temperature to the energy of the parton, 
v T/E, the contribution coming from quark-quark scattering can be reduced 
to [113] 



fdW_\ Qq 
\dx) 



2(AixfN f a 2 s f d 3 k / c (k) f d*k> l-ft(U 



hard 



3v 

x G(q-q*) 



r d A k £(k) r d 6 i 

J j2^~k~J (2^ 



UJ 



(uj 2 — q 2 Y 



2(k- v-k) 2 + 



-S(uj 



(- 2 - <f) 



q) 



, (7.22) 



after performing the Dirac traces and evaluating the sum over spins while the 
contribution coming from quark-gluon scattering gives 



aw 

dx 



(47r)X f d 3 k n ? (k) f d 3 k' l+n f (k' 



hard 



2v 



X UJ 



(2tt) 3 

(1-v 2 ) 2 
(k- v-k) 2 



k J (2tt) 3 k' 

(k-vk) 2 + ±=f(uj 2 



+ 8 



5(cu-vq)e(q-q*) 



(uj 2 - q 2 )' 



(7.23) 



Here f^(k) and n^(k) are the anisotropic versions of the tree-level Fermi-Dirac 
and Bose-Einstein distribution functions at zero chemical potential and uj = k' — k 
while q = k'—k. Note also that q* acts as IR cutoff for the q integration. Since the 
integrand is odd under the interchange k <-> k' the terms involving the products 
/^(k)/^(k') and n^(k)n^(k') vanish since they are symmetric. Redefining the 
origin of the k' integration so that it becomes an integration over q one obtains 

'd^\ Q9 



dx / 



hard 



16N f aj r d A k /g(k) 
~3v J (2tt) 3 k 
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UJ 
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4a 2 f d 3 k n ? (k) 
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5(uj 
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k) 2 + 



1-v 2 



[UJ^ 
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(7.24) 
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Figure 7.4: Tree-level Feynman diagrams for the scattering processes Qq 
(first diagram) and Qg — > Qg (remaining diagrams). 



Qq 



where now uj = |q + k| — k. 

Choosing v to be the z-axis for the q and k integration one can rewrite the 
delta-function as 



S(\q + k| — k — v • q) 



5(0 9 -0 o )e(A; + v-q)2|q + k| 



q\J Ak 2 sin 2 9k sin 2 9 q — (g(l — v 2 cos 2 9 q ) + 2 cos 6 q (k cos 9~k — kv)) 2 
where 0o is the solution of the equation 

q(l — v 2 cos 2 9 q ) + 2 cos9 q (k cos 9k — kv) 



,(7.25) 



COS! 0o - <Pk 



2k sin 9k sin 9 g 



(7.26) 



Including a factor of 2 because of the symmetry (p ^ 2tt — O the g integration 
is straightforward. Moreover, the q integration can be done by scaling k = qz as 
well as uj = qvcos9 q , since then 



n € (k) -> y/T+£n(qzJ 1 + £(k • n) 2 ) 



(7.27) 



and one only needs to consider the integrals over Fermi-Dirac and Bose-Einstein 
distributions, 

dqqf(q) = ^ + ?W (-*)) - U 2 (- e W (-x)) = ^ 



x 2 T 2 



xT 
oo 

xT 



X, 



dq q n(q) 



-x\n(l — exp (—x)) + Li 2 (exp (—x)) 
x 2 



Fo(x) 



where x = %^-\/\ + £(n x sin# fc cos0fc + n z cos9 k ) 2 , with n 2 = nv and 1 
The hard contributions to the energy loss then take the form 



(7.29) 
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Figure 7.5: Hard contribution to the isotropic energy loss for Nf = 2 (full line) 
as a function of q* for v = 0.5 compared to the Braaten-Thoma result (dashed 
line). 
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(7.30) 



where T denotes the unwieldy expression 



T = Q(z+v cos9 q 



Q(Az 2 sin 2 9 k sin 2 9 q - (1 - v 2 cos 2 9 q + 2 cos 6»^(cos 6» fc - v)) 2 ) 



y / 4^ 2 sin 2 fc sin 2 ^ — (1 — f 2 cos 2 6 q + 2 cos^ g 2;(cos ^ — w)) 2 

(7-31) 

The remaining integrations have to be performed numerically. 



7.2.1 Behavior of the hard part 

Similar to what has been done for the soft part it is possible to gain some in- 
sight on the parametric dependence of the hard contribution by considering the 
anisotropy to be weak (£ 1) [114]. One can then compare the unexpanded 
isotropic (£ = 0) result to the Braaten-Thoma result, shown in Fig. 7.5. As one 
can see, the two results are identical for small q* = q*/T while for larger q* the 
Braaten-Thoma result becomes negative and the unexpanded result Eq.(7.30) is 
positive for all values of q* . 



7.3. ISOTROPIC RESULTS 
a 



b 
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Figure 7.6: Isotropic energy loss of a beauty (a) and charm quark (b) as a function 
of momentum p for T = 250MeV and a s = 0.3. Shown are the respective results 
from Ref. [117] (dotted line) and Ref. [113] (dashed line) as well as the evaluation 
of Eq.(7.32) (full line) with a variation of q* (gray band). 

7.3 Isotropic results 

In the isotropic limit, the total collisional energy loss of a heavy parton is obtained 
by adding Eq.(7.21) and Eq.(7.30) with f = 0, 



dW 

dx 



dW 

dx 



+ 



soft,iso 



dW 

dx 



hard,iso 



(7.32) 



It is a function of the strong coupling a s , the particle velocity v, the temperature 
T and the momentum separation scale q* . As has already been discussed above, 
the q* dependence of the result is found to become weak for small values of the 
coupling limit a s (similar to what has been found in QED [114]), corresponding 
to the original result by Braaten and Thoma [113]; for larger values of the 
coupling one fixes q* = q pms using the principle of minimum sensitivity, 



d 
dq* 



dW 

dx 



= 0. 



(7.33) 



<?*=<?; 



Note that — (-7^). , always serves as a lower bound on the result for the 

energy loss. To get an estimate of how strongly the result depends on this special 
value of q*, one can e.g. vary q* ms by a certain factor c q * and evaluate the energy 
loss at <?p ms c 9 * and (?p ms /c ? *. This factor c q * is in principle arbitrary, but should 
be such that the resulting q* is neither much smaller than mo, nor much bigger 
than 2nT; in the following, I have chosen c q * = 2 to be consistent with what has 
been done in QED [114]. 

In Fig. 7.6a,b various results for the energy loss of a heavy parton in an 
isotropic quark-gluon plasma are compared: shown are the results from BJORKEN 
[117], Braaten and Thoma [113], as well as Eq.(7.32) at q* = g* ms together with 
its variation using c q * = 2. For the compilations a temperature of T = 250MeV 
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and a coupling constant of a s = 0.3 was adopted; Fig. 7.6a shows the energy loss 
of a beauty quark with mass Mq ~ 5GeV while in Fig. 7.6b the energy loss of 
a charm quark with mass Mq ~ 1.5GeV is plotted, both as a function of their 
momenta 

" = 7r=^ (7 ' 34) 



7.3.1 Limitations 

Since Eq.(7.32) has been derived for an infinitely heavy parton, it breaks down 
for thermal velocities v ~ y/T '/Mq because a quark with v — can only gain 
energy from collisions with other particles in the heat bath. A semi-quantitative 
estimate of the velocity where the energy loss becomes negative has been done 
in Ref. [113] by repeating the above calculation in the limit v — > and for 
weak couplings, finding v ~ a/3T/Mq (which corresponds to p = 1.5GeV and 
p = 2.1GeV for charm and beauty quarks, respectively). Similarly, Eq.(7.32) also 
breaks down for ultrarelativistic energies E ^> Mq/T, with the cross-over energy 
having been determined to be -E CT oss — 1.8Mq/T (corresponding to v > 0.995 
for both charm and beauty quarks). However, note that the reason that the 
Braaten-Thoma result turns negative for a momentum of p ~ 5.7GeV for the 
beauty and p ~ 1.71GeV for the charm quark (as is indicated in Fig. 7.6a,b) is 
not due to this physical reason but rather due to a failure of the extrapolation 
from the weak coupling limit to realistic couplings [113]. For the unexpanded 
result Eq.(7.32), this unphysical behavior does not occur and one can therefore 
expect it to be valid for velocities down to the original estimate v ~ a/3T/Mq. 

Finally, it should be noted that no estimate on the next-to-leading order 
(NLO) corrections to the energy loss has been made here. Therefore, it should be 
kept in mind that for QCD with large realistic coupling the inclusion of these NLO 
corrections might give energy loss results that are not covered by the variations 
of q* in the leading-order result, so these variations should be interpreted with 
care. 



7.4 Anisotropic results 

The full result for the heavy parton energy loss in an anisotropic system with 
strength £ is given by adding the soft contribution from Eq.(7.18) and the hard 
contributions from Eq.(7.30), 



dW 

dx 



dW 
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+ 



soft 



dW 

dx 



hard. 



(7.35) 



It depends on q*,T and v as was the case for the isotropic energy loss, but in 
general also on the angle of the particle direction with respect to the anisotropy 
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Figure 7.7: Anisotropic energy loss 
at v = 0.5, £ = 1.0 for a s = 
0.2, 0.3, 0.4 (dotted, full and dashed 
line, respectively) as a function 
of the angle 9 (normalized to the 
isotropic result). 
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Figure 7.8: Anisotropic energy loss 
for f = 1.0, a s = 0.3 and 9 = (full 
line) compared to 9 = tt/2 (dashed 
line) as a function of velocity (nor- 
malized to the isotropic result). 



vector cos 9 = v • n (sketched in Fig. 7.1). The angular dependence of the energy 
loss at £ = 1 normalized to the isotropic energy loss Eq.(7.32) is shown in Fig. 7.7 
for v = 0.5 and different couplings a s = 0.2, 0.3, 0.4; as can be seen, for all cou- 
plings considered the energy loss turns out to be smaller along n than transverse 
to it. However, as already noted in the earlier discussion on the soft contribution 
to the energy loss, this dependence may change as a function of the velocity v, 
as can be seen in Fig. 7.8, where the energy loss Eq.(7.35) evaluated at 9 = 
and 9 = n/2 has been plotted for a s = 0.3. From this figure it becomes clear 
that the energy loss is peaked at 9 = ir/2 for velocities smaller than vo while it is 
peaked at 9 = for v > vq, where t>o — 0.7 for £ = 1 and a s = 0.3 (see also [118]). 
Note that in Figs. 7.7, 7.8 the variational bands correspond largely to those of 
the isotropic results (about 40% higher). 

For simplicity here it was assumed that quark and gluon distribution func- 
tions have momentum-space anisotropies of the same strength. In general, this 
assumption is quite probably not fulfilled, so a more realistic result for the energy 
loss will depend on the quark and gluon anisotropies separately. 



7.5 Effects of finite chemical potential 

Taking finite chemical potential into account, the soft contribution Eq.(7.18) 
remains unchanged in form with the only difference being the change of the 
Debye mass 

ml 4 TO< ( 6 -±flr- + !g) , ,7.36, 
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in the structure functions according to Eq.(2.4). For the hard contribution, only 
the quark-quark scattering term changes, so one has to replace 



in Eq.(7.22). Again, all terms involving products of distribution functions can 
be shown not to contribute to the energy loss, so one recovers Eq.(7.30) with the 
function Fi(x) being replaced by 



which can be calculated analytically similar to Eq.(7.28). A subsequent evalua- 
tion of the collisional energy loss at finite chemical potential in the weak coupling 
limit has been done in Ref. [119]. There it was found that the effects played by 
the quark chemical potential are rather small unless \ijT gets large. 

7.6 Summary 

In this chapter I have calculated the complete leading-order collisional energy 
loss of a heavy parton propagating through an anisotropic quark-gluon plasma. 
The calculation was based on results for the collective modes in an anisotropic 
plasma of chapter 6 and it was demonstrated that the unstable modes found for 
such a system would in principle cause the result to be divergent. However, it 
has been shown that the collisional energy loss is protected by the mechanism of 
dynamical shielding, for which a proof in two analytically tractable regimes (weak 
and extremely strong anisotropy) was given. As a side result, in the isotropic limit 
the original result for the energy loss was shown to have a correction which cured 
the problem of unphysical energy loss results for heavy partons. 

The results were applied for beauty and charm quarks and the applicability of 
the calculation was discussed. Also, it was shown that for anisotropic systems the 
collisional energy loss has an angular dependence which is expected to increase for 
larger couplings and stronger anisotropies, possibly having phenomenological im- 
plications. Finally, a refinement of the calculation involving different momentum- 
space anisotropies for the quark and gluon distribution functions was suggested 
and the effects of finite chemical potential were discussed. 



/ € (k)(i - / fi (k0) - + f-tto) (i - (iW(k) + f-*M) (7-37) 




(7.38) 



Chapter 8 

Conclusions and outlook 



In this work, I have investigated the quasiparticle description of certain observ- 
ables in the hot and dense quark-gluon plasma. More precisely, the quasiparticle 
excitations and in general the collective modes of systems with isotropic as well 
as anisotropic momentum-space distribution functions were analyzed in the HTL 
approximation, with a special emphasis on the gluonic excitations. 

For isotropic systems, an HTL quasiparticle model for the thermodynamic 
pressure based on the two-loop ^-derivable approximation for the entropy was 
proposed, which improves on an existing simpler model by taking into account 
the full momentum dependence of the HTL self-energies. It was shown that these 
quasiparticle models can be used to accurately describe two-flavor lattice results 
for the pressure, entropy and quark-number susceptibilities at zero chemical po- 
tential by adopting an effective running coupling including two fit parameters. 
An extension of these quasiparticle models to finite chemical potential (naturally 
given by a flow equation for the strong coupling as a consequence of the sta- 
tionarity condition for the quark and gluon propagators) was also shown to be 
consistent with independent lattice studies in the region of applicability of the 
latter. Therefore, these results demonstrate that quasiparticle models can qual- 
itatively reproduce lattice results for the major bulk thermodynamic quantities 
above the phase transition, at least when supplied by a phenomenologically in- 
spired effective strong coupling as an input. Quantitatively the differences with 
respect to lattice studies were found to be small but non-vanishing, although it 
is probably fair to say that currently the variations in the latter due to different 
implementations of fermions and incomplete continuum-extrapolations are com- 
parable in magnitude. Encouragingly, the inclusion of the plasmon effect which 
has devastating effects on a strictly perturbative calculation for the pressure turns 
out to lead to small corrections when incorporated through a quasiparticle model, 
as has been shown here. 

Therefore, quasiparticle models in general and the HTL model in particu- 
lar provide a simple, physically intuitive and rather accurate description of the 
quark-gluon plasma, both at zero and also at non-vanishing chemical potential. 
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Moreover, motivated by this success at small chemical potential, estimates of 
quantities at large chemical potential, such as the critical chemical potential at 
very small temperature, can and have been calculated. Although these results 
represent only estimates based on the physics at small \i and are therefore likely 
to change somewhat once techniques which are more accurate for the high den- 
sity regime are developed, the quasiparticle model calculations give precise and 
robust predictions in contrast to lattice studies which (as yet) cannot reach this 
domain of high density at all or strictly perturbative results which suffer from 
large renormalization scheme dependences. 

However, further refinements are possible: as a first step, a complete determi- 
nation of the next-to-leading order corrections to the self-energies would eliminate 
the dependence of the NLA model on the free parameter c\, allowing a precise 
inclusion of the full plasmon term in the model, thereby greatly increasing its 
accuracy. The next (ambitious) step would be to aim for a fully self-consistent 
calculation of the one-loop self-energies, which in turn would make the quasipar- 
ticle model completely self-contained, since both the form of the gap-equations as 
well as the running coupling would result from such a calculation. Unfortunately, 
while the first step is fairly straightforward, the second step probably still has 
to await progress concerning gauge-dependencies and the renormalization pro- 
cedure of e.g. the two-loop ^-derivable approximation for QCD. Nevertheless, 
the results of such calculations would be highly interesting since it is possible 
that knowledge about such fully non-perturbative quasiparticle-like excitations 
might be sufficient to offer something that is very close to a full description of 
the isotropic quark-gluon plasma. 

For anisotropic systems, however, the situation is fairly different: the reason 
for this is that systems with an anisotropy in momentum space possess unstable 
modes that correspond to exponentially growing gauge field amplitudes. The 
presence of these unstable modes is believed to have important consequences 
both on the dynamical evolution of the quark-gluon plasma as well as on more 
practical issues such as the de-facto breakdown of calculations using a pertur- 
bative framework, as has been demonstrated. However, owing to the fact that 
in QCD the (re-) discovery of these instabilities has only occurred very recently, 
many important questions remain largely unanswered and indeed are subject to 
active ongoing research. For example, there are indications that certain observ- 
ables exist that are protected from the breakdown of perturbative calculations, 
as has been shown explicitly here for the collisional energy loss of a heavy par- 
ton. However, it is currently unknown if there are any other observables that are 
either dynamically shielded or otherwise protected from the singularities coming 
from the unstable modes, so further studies in this direction are mandatory to 
have a better understanding of the effects of the QCD plasma instabilities. Also 
the saturation mechanism of the instabilities deserves closer investigation: while 
for QED non-perturbative effects on the hard particles are responsible for halting 
the growth of the unstable modes, for QCD the non-Abelian interaction between 
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the soft fields might saturate this growth already earlier. Moreover, a detailed 
understanding of the role that the instabilities play in modifying the bottom- 
up thermalization scheme is still missing and would probably be very important 
for the early system evolution following a heavy-ion collision. Summarizing, the 
full implications of unstable modes in the anisotropic quark-gluon plasma have 
not yet been worked out, but current knowledge suggests that these might have 
(measurable?) consequences on the physics tested at RHIC and especially the 
LHC once it becomes operable, since there the momentum-space anisotropies are 
probably very strong. 

In conclusion, calculations based on the analysis of collective modes offer 
qualitative as well as quantitative predictions relevant for QCD at high temper- 
ature (and non-vanishing chemical potential) both for isotropic and anisotropic 
systems. However, whereas for isotropic systems quasiparticle models for thermo- 
dynamic quantities are well-developed and have been tested in many occasions, 
anisotropic systems are less well understood and their theoretical treatment has 
only been tested in the case of QED. Moreover, current calculations (although 
non-perturbative in character because of implicit resummations) are essentially 
limited to including only leading-order correction terms completely, so the quan- 
titative (as well as qualitative?) effects of higher order corrections might still 
be substantial for QCD where the coupling gets large. Nevertheless - pending 
the development of new methods - quasiparticle descriptions of hot and dense 
QCD allow predictions in a range of situations where other approaches from first 
principles become either ambiguous, very complicated or break down altogether, 
thus making the former an invaluable tool in the study of the new form of matter 
called the quark-gluon plasma. 
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Appendix A 



Table of symbols 



For convenience, I include here a short explanation of some of the symbols used 
in the main text (roughly in order of appearance): 



Symbol 


Meaning 


T 


Temperature 


fi 


Quark chemical potential 


T 

A c 


Critical temperature at /i — 


p 


Thermodynamic pressure 


Po 


Stefan-Boltzmann pressure 


Ho 


Chemical potential at T ~ where p = 


He 


Chemical potential at T ~ where p = p{T — T c ,/i — 0) 


p, 


Renormalization point in MS scheme 


A MS 


Renormalization scale in MS scheme 


s, n, e 


Entropy, quark number and energy density 




Quark number susceptibilities 


a s 


Strong coupling constant 


N f 


Number of quark flavors 


m D 


Debye mass 


T s , A, B 


Quasiparticle model fit parameters 


z 


Anisotropy parameter (see chapter 6) 




Anisotropic self-energy structure functions 


Q 


Momentum separation scale 


dW/dx 


Collisional energy loss 
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Appendix B 

Coefficients for the HTL coupling 
flow equation 



In order to rewrite Eq.(3.17) into the form of Eq.(2.26), quite some algebra is 
necessary. For convenience, I therefore reproduce here the final coefficients a^, a M 
and b which may be directly used to solve the flow equation. I find 

2NN f C 



where 
A = 



B 



ax 



Ana 



s,eS 



N 2 - IT' 
_A_ _ 2NN f B 

~ jV 2 - IT 2 ' 
2NN f C 47ra s , cff N 2 - 1 
N 2 — IT M 2 8N~ 
2NN f B 4vra SiCfr N 2 — 1 /j, 



_ A ATca S!eS N f 

rrl 2 2 f 1 

1 z m D 7T Z 



d 3 k 

(2^)3 

■ k 



N 2 -IT 2 M 2 8N n 2 ' 

2(k 2 - LU 2 )(lmU T ) 3 



did dn(uj)(—uj) 



2 f 

* Jo 

■2E T dn(E T ) z?2 
d 3 k 



[(k 2 -oj 2 + Ren T ) 2 + (Imn T ) 2 ] 2 
A; 2 (Imn L ) 3 



[{k 2 + Ren L ) 2 + (Imn L ) 2 ] 2 
E 2 L dn{E L ) 



E 2 -k 2 



El - k 2 



m 2 D E 2 — (E T — k 2 ) 2 

k 



\k 2 -E 2 +m 2 D \ 
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n Jo 



x 
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+ 
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[(A; + cj + Re£_) 2 + (ImE_) 2 ]' 



89 



90APPENDIX B. COEFFICIENTS FOR THE HTL COUPLING FLOW EQUATION 



+ (E + - k) [df + {E + ){-E + + p) - df_(E + )(E + + /i)] 
+ {E- + k) [df + (E_)(-E_ df_(E_)(E_ + /i)] 



E\ - k 2 

2M 2 
E 2 _-k 2 

2M 2 



(B.3) 



C = 



d 3 k 



7T 



- / dw (-df + (u) + df.(u)) 



(k — w)(ImS_ 



[(k — u + ReX+) 2 + (ImS + ) i 
(A; + cj)(ImS_) 3 



+ - / duo [-df + (u) + 

Wo [(A; + cc> + ReX_) 2 + (Im£_) 2 ] 2 

i-df + (E + ) + df_(E + )} 
+ (E^ + k) [-df + (E_) + df4E_)]}. 

In the above equations I have used the abbreviations 

— expu/T 



dn(ui) 
df ± (uj) 



(expw/T- l) 2 

-exp {u±n)/T 
(exp(cu±/x)/T+ 1) 



2 ' 



(B.4) 



(B.5) 



Furthermore, E T ,E L and E + ,E_ are the dispersion relations E T (k), E L (k) and 
E + (k), E_(k) for the bosonic and fermion quasiparticles shown in chapter 2, and 
I1 T , n^, E + , £_ are the bosonic and fermionic HTL self-energies also presented in 
this chapter. 



Appendix C 

Analytic expressions for structure 
functions 



In this appendix I collect the integral and analytic expressions for the structure 
functions a, j3, 7, and 5 defined in Eq.(6.31). Choosing n = z and k to lie in the 
x-z plane (k x / k z = tan 9 n ) one has vn = cos 9 and v k = v k x cos (f> sin 9+vk z cos 9. 
Using this parameterization the integration for all four structure functions 
defined by the contractions in Eq.(6.31) can be performed analytically. 



m 



D 



k 2 n 2 
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uj — k, cos 9 
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uk z . 
k 2 ■ 
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m 2 D uj 2 
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2tt a - 6cos0 + ie \/a + 6 + ieVa -b + ie 
When a and b are real-valued i? can be simplified to 

sgn(a)0(a 2 - b 2 ) iQ{b 2 - a 2 ) 
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(C.3) 
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with Q(x) being the usual step-function. Note that the remaining integration 
over 9 can also be done analytically but the results are rather unwieldy so I do 
not list them here. 

Static Limit 

In the limit ui — > it is possible to obtain analytic expressions for all four structure 
functions. The results for m a and mp defined in Eq.(6.38) are 





arctan 



,(C7) 




2 ( v / e + (l+0arctan v ^)(k 2 + ek2) 




(C.8) 



with similar results for mi and m\. 
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